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Overview: 


During  the  last  contract  period  we  have  published  two  papers,  we  have 
one  "in  press,"  with  a  third  ready  for  submittal  and  a  fourth  in  the 
late-draft  stages  (see  enclosures).  The  research  culminating  in  these 
papers  has  been  supported  entirely  or  in  part  by  ONR  contract 
N0001A-75-C-0539.  We  also  have  developed,  modified,  applied  and 
documented  a  three-axis  image-analysis  system  for  analyzing  holograms  of 
dynamic  particles.  A  description  of  the  image-analysis  hardware, 
algorithms  for  particle  sizing  and  the  accuracies  attained  can  be  found  in 
Payne,  Carder,  and  Steward  (1983;  see  attachment).  A  description  of  the 
nine  holographic  image  analysis  subroutines  for  controlling  a  three-axis, 
2-mlcron  resolution  positioning  system,  particle  sizing,  and  performing 
data  management  functions  is  also  provided  along  with  program  listings. 

Synopsis  of  Results 

The  two  papers  published  last  summer  resulted  from  earlier  N'avy 
contracted  work.  They  are: 

Fanning,  K.  A.,  K.  L.  Carder  and  P.  R.  Betzer,  1982.  Sediment 

resuspension  by  coastal  waters:  a  potential  mechanism  for 
nutrient  re-cycling  on  the  ocean's  margins.  Deep  Sea  Research 
mSA):  953-965. 

Carder,  K.  L.,  R.  G.  Steward  and  P.  R.  Betzer,  1982.  In  situ 

holographic  measurem.ents  of  the  sizes  and  settling  rates  of 
oceanic  particulates,  Geophysical  Research  87 (C8) : 

5681-5685. 

In  the  first  paper  we  demonstrated  that  very  large  light-attenuation 
coefficients  (e.g.,  C^  >3.3  m  S  measured  near-bottom  (■^30m)  during  a 
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winter  frontal  passage  over  our  three-day  time  series  station  south  of 
Mobile  Bay  were  accompanied  by  nutrient  enrichments  of  as  much  as  10-fold 
over  concentrations  found  in  the  overlying  waters  or  between  storms.  We 
provided  evidence  that  sediment  resuspension  may  significantly  accelerate 
nutrient  recycling  on  continental  margins. 

The  second  paper  documents  the  first  ^  situ  oceanic  measurements  of 
individual  particle  settling  velocity  parameters.  They  were  obtained 
using  a  holographic  technique  in  a  vertically  damped  sediment  trap/float 
array  in  the  western  Sargasso  Sea.  To  our  surprise  there  were  significant 
numbers  of  large  (15-45  ym)  heavy  minerals  of  Saharan  dust  origin.  Such 
particles  have  been  heretofore  never  reported  in  the  aerosol  literature  so 
far  (>5000  km)  from  the  Sahara  Desert.  The  holographically  determined 
sizes,  shapes,  and  Stokes-derived  densities  matched  those  determined  by 
scanning  electron  microscopy  analysis  (SEM/EDXA)  and  particle  compositions 
of  samples  collected  in  the  sediment  trap. 

This  research  activity  also  led  to  three  papers  which  are  in  various 
stages  of  publication; 

Doyle,  L.  J.,  K.  L.  Carder,  and  R.  G.  Steward,  The  hydraulic 

equivalence  of  mica.  Sedimentary  Petrology,  (in  press). 

Payne,  P. ,  K.  L.  Carder  and  R.  G.  Steward,  A  system  to  catalog  and 

analyze  holographic  images  of  dynamic  ocean  particles,  final 
draft  to  be  submitted  to  Applied  Optics . 

Carder,  K.  L.,  R.  G.  Steward,  P.  R.  Betzer,  D.  Johnson  and 

J.  Prospero,  Chronology  of  an  aeolian  input  event  to  the 
Sargasso  Sea,  to  be  submitted  to  Geophysical  Research. 

The  first  paper  reports  the  findings  of  a  laboratory  experiment  using 
the  holographic  nicroveloclmeter .  The  purpose  of  this  experiment  was  to 


2 


determine  the  settling  dynamics  of  sand-  to  silt-sized  mica  flakes.  Mica 
has  been  classically  considered  to  have  a  settling  velocity  relatively 
lower  than  other  minerals  of  equivalent  size  and  density.  This  was  based 
upon  the  association  of  sand-sized  mica  and  clay-sized  minerals  in  various 
environments.  Until  our  experiment,  there  had  been  no  quantitative 
evidence  to  support  this  correlation.  The  holographic  microvelocimeter 
provided  the  means  of  accurately  viewing  and  measuring  the  settling 
behavior  and  velocities  of  individual  mica  flakes.  These  data  established 
that  mica  is  the  hydraulic  equivalent  of  quartz  spheres  between  4  to  12 
times  smaller  in  diameter. 

This  experiment  was  conducted  using  dual  beam  holography;  i.e. 
simultaneous  holograms  were  collected  along  the  vertical  and  horizontal 
axes.  These  data  showed  that  the  mica  flakes  seldom  settled  in  either  the 
broadwise  or  edgewise  orientation.  Also  the  flakes  usually  maintained  a 
stable  orientation  and  rarely  rotated  further.  The  data  allowed 
comparison  between  two  different  formulations  of  Stokes  equations  modified 
for  shape  and  settling  orientation.  This  comparison  showed  that  for 
flakes  where  the  diameter  is  approximately  50  times  the  thickness  there 
was  little  difference  between  the  two  sets  of  shape  modified  equations. 
For  thicker  particles,  the  orientation  factors  in  one  set  of  equations 
allowed  better  accuracy  in  predicting  the  settling  velocity. 

The  second  paper  (see  Appendix)  describes  the  image  analysis 
hardware,  algorithm  concepts,  and  algorithm  accuracies  for  a 
micro-computer  controlled  hologram  positioning  and  particle  analysis 
system.  Particle  size  (area,  perimeter),  shape,  settling  velocities,  and 
Stokes  densities  are  measured/calculated  with  the  system  with  accuracies 
of  better  than  14%  for  high-contrast  particle  images  and  better  than  ±10% 
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for  low-contrast  holographic  particle  images  with  optical  noise  (speckle 
and  particle  diffraction  patterns)  superimposed.  Techniques  for  improving 
signal  extraction  from  the  noisy  images  are  also  discussed. 

The  third  paper  discusses  development  of  a  three-week  time  series  of 
Sahara-derived  aerosol  concentrations  from  the  particle  data  collected  by 
the  holographic  particle  velocimeter  trap  data.  In  essence,  a  particle 
settling  model  (particle  densities  and  shape  parameters)  consisting  of 
eight  particle  classes  (e.g.,  quartz,  hematite,  rutile,  etc.),  was 
calibrated  by  fitting  the  aerosol  concentration  data  derived  from  taking 
the  particles  from  the  sediment  trap  (30  m  depth)  back  in  time,  to  aerosol 
concentration  data  collected  by  Joe  Prospero  on  Virginia  Key.  This  in 
effect  was  a  natural  sedimentation  experiment  with  a  30-meter  "settling 
tube."  It  produced  Stokes  particle  densities  for  each  mineral  class. 
These  matched  almost  precisely  our  holographically-derived  densities 
(Carder,  e^  ,  1982)  for  heavy  minerals  (e.g.,  5.10  vs.  5.20  for 
"hematite,"  4.10  and  4.16  vs.  4.20  for  "goethite"  and  "rutile"),  and  for 
non-heavy  minerals  (e.g.,  from  2.55  to  2.72  versus  2.19  to  2.99).  Since 
the  non-heavy  minerals  included  clays,  quartz,  and  Iron/titanium 
containing  silicates,  the  variance  of  density  values  for  these  non-heavy 
mineral  classes  would  be  expected  to  be  larger  than  for  the  dominant  heavy 
mineral  classes. 

Thus  success  of  this  particle  dynamics  model  for  aeolian  input  will 
permit  us  to  now  look  at  fluxes  of  this  component  as  a  tracer  through  the 
ocean  system  (e.g.,  removal  and  repackaging  by  filter  feeders)  in  future 
papers.  With  data  collected  over  5  days  trom  traps  at  30  m,  100  m  and 
300  m,  development  of  aerosol  time  series  for  the  Sargasso  of  up  to  2 
months  should  be  possible.  This  would  provide  us  the  capability  of 
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gathering  time  series  data  on  aerosol  inputs  to  the  ocean  in  remote 
settings  without  the  expense  of  maintaining  a  ship  on  station  for  long 
periods  of  time.  The  investigation  of  the  temporal  and  spatial  coherence 
of  aeolian  Inputs  to  the  sea  might  be  useful  future  outgrowths  of  this 
line  of  research  as  well  as  to  provide  increased  understanding  of  the  real 
and  "optical  residence  times"  that  Saharan  and  Gobi  Desert  dust  have  with 
respect  to  light  propagation  through  various  regions  of  the  ocean.  Our 
Sargasso  Sea  data  indicate  chat  Saharan  dust  particles  smaller  than  about 
5.0  Urn  diameter  have  a  residence  time  in  the  upper  30  m  of  more  than  two 
weeks  (8  weeks  for  2.5  Urn  particles).  Since  during  the  summer,  Saharan 
dust  input  events  occur  every  7  to  14  days,  an  optical  signal  may  persist 
in  the  upper  layers  for  days  after  an  aerosol  cloud  has  passed  by. 
Whether  or  not  this  aeolian  residue  has  a  significant  effect  on  interpre¬ 
tation  of  remotely  sensed  optical  data  or  optical  signal  propagation 
awaits  future  study. 

Other  rescrlts  apparent  from  the  data  of  Carder  a^.  (in 

preparation)  are  as  follows: 

1.  Aerosol  size  distributions  and  mineral  types  change 
significantly  v/ith  the  phase  of  the  aerosol  clouds: 

a.  A  disproportionately  larger  fraction  of  the  particles  are 
heavy  minerals  and  "giant"  particles  (  10  Pm  diameter)  in 
the  leading  edge  of  the  cloud  than  at  the  trailing  edge. 
Aeolian  particles  in  excess  of  60  urn  diameter  were 
observed. 

b.  Conversely,  the  trailing  edge  has  a  paucity  of  large,  heavy 
particles. 
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c.  A  far  larger  fraction  in  sununer  dust  clouds  are  heavy 
mineral  particles  than  in  winter,  suggesting  that  the 
summer  source,  which  is  considerably  north  of  the  winter 
harmatan  source,  is  a  more  titanium/hematite/goethite-rich 
region. 

IMAGE  ANALYSIS  SOFTWARE 

The  image  analysis  software  developed  in  collaboration  with  Aztec 
Computer  Engineering  for  this  contract  is  composed  of  nine  subroutines. 
The  source  code  of  each  is  listed  in  the  Appendix.  With  this  software  the 
analysis  of  the  holographic  image  is  now  semi-automated.  The  computer 
permits  precise  registration  of  the  holograms,  high  resolution  scanning  of 
the  video  images  to  measure  the  particle  area  and  perimeter,  calculation 
of  the  particle  fall  distance  between  frames  for  velocity,  and  data  base 
management.  These  software  modules  are  described  below. 

1.  MAIN  The  main  routine  controls  the  flow  of  the  other  modules 
through  a  menu,  calculates  particle  settling  velocity,  and  writes  the  data 
onto  floppy  disks, 

2.  HELP  This  is  a  help  routine  which  can  be  entered  at  various 
places  in  the  main  routine  to  explain  the  program  operation,  operator 
inputs,  and  computer  responses. 

3.  FMOVE  Frame  move  is  a  routine  to  control  the  hologram 
positioning  in  three  dimensions  using  the  computer  joystick  and/or  the 
keyboard. 

k.  STAGE  This  is  the  interface  subroutine  that  interprets  the 

output  of  the  joystick  and/or  keyboard  and  moves  the  XYZ  stage  motors. 
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5.  JSINP  Joystick  input  is  an  interfacing  subroutine  which  outputs 
a  digital  number  from  the  computer  joystick  to  be  used  in  control  of  the 
XYZ  stages  or  the  video  cursor. 

6.  PSCAN  Particle  scan  is  a  routine  that  measures  the  area  and 
perimeter  of  a  video  image.  This  algorithm  operates  by  looking  for  video 
pixels  of  higher  intensity  than  an  operator  selected  threshold. 

7.  TVSCN  Television  scan  is  a  subroutine  of  PSCAN  which  allows  the 
operator  to  move  the  video  cursor  to  the  center  of  the  video  screen, 
particle  center,  and  particle  extremes.  After  setting  the  extremes  of  the 
particle,  this  routine  sums  up  the  area  by  counting  all  pixels  greater 
than  the  selected  threshold  and  sums  up  the  perimeter  from  the  number  of 
the  rows  and  columns  of  the  scanned  area. 

8.  TRACK  This  is  the  particle  edge  tracking  routine  which 
automatically  tracks  around  the  edge  of  the  video  image.  This  algorithm 
tests  for  an  operator-selected  video  amplitude  gradient  that  approximates 
the  particle  edge  conditions. 

9.  TVTRK  The  television  tracking  routine  is  a  subroutine  of  TRACK 
which  moves  the  video  cursor  around  a  circle  at  each  XY  position  on  the 
particle  edge.  In  moving  around  this  circle,  the  video  amplitude  is 
tested  to  see  if  two  adjacent  perimeter  points  on  the  circle  have 
difference  in  video  amplitude  which  exceeds  the  selected  gradient.  If  a 
complete  circle  is  made,  or  if  the  gradient  is  exceeded  then  the  routine 
returns  the  XY  coordinates  of  the  perimeter  point  to  be  used  as  the  center 
of  the  next  scan. 
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MAIN  PROGRAM  FOR  THE  AZTEC  HOLOGRAPHIC  IMAGE  ANALYSIS  SYSTEM 
REWSION  0.30,  5  SEP  1982 

LOGICAL  ESC. AST, ANS,CHG,FILE(  U  ),HA«Ei 3 1.DATE'  3 l.UORIiSl  44 1 
INTEGER  RMAX,FMAX.FRAME,PART,FR1.FR2 
INTEGER  FRM(  10).PRT(  10),PX(  10).PY<  10).PZ(  10).3EX(  10),EEi(  10) 
INTEGER  SAREA(10),SPERIM(IO),SHID£(10).STALL;10) 

INTEGER  TAR£A(  10 ),  TPERIHl  10 ),  TLEN!  10 ),  TX(  10.100 ),  TY<  10. 100 ) 
INT£G£R  HLP.FX,FY,FZ,SX,SY,SA,SP,SH,SH,TA,TP,TL,T«(  lOOi.YHI  100 
DIM£NSION  TIMeI 10 )  ’  ’ 

EXTERNAL  HELP. FHOVE.PSCAN. TRACK 
tXTERNAL  JSINP, STAGE, TVSCN. TVTRK 

COMMON  ESC,AST,HLP,FX,FY,F2,SX,SY,SA,3P,SW,SH,TA,TP.TL.7H.TH 

DISPLAY  TITLE  AND  OVERVIEW  SCREEN  <HLP=1) 

£3C=:7 
AST =42 
HLP=1 
CALL  HELP 


C  INITIALIZE  SUBROUTINE  PARAMETERS 

FX=0 

FY=0 

FZ=0 

?X=128 

SY=123 

3A=0 

3F'=0 

SU=0 

SH=0 

TA=0 

TP=0 

41.100 

th(:j)=o 

105  CONTINUE 

c 

C  SELECT  TASK  OR  DISPLAY  H.L.T.X  SCREEN  (HLP=2) 

C 

110  HRITEvUllOO) 

1100  FORMAT!'  SELECT  (L)OAD  OR  (T)EST:  '' 

REAIKl.mOANS 
1110  FORMAT!  lAl) 

HLP=2 

IF!ANS.EQ.72)CALL  HELP 
IF!ANS.E0.76)G0T0  200 
IF1ANS.E0.84)G0T0  120 
IF!ANS.Ei}.3e)G0T0  190 
GOTO  110 
C 

c  SELECT  TASK  OR  DISPLAY  H,F,P,T,X  SCREEN  'HLP-S) 

i:0  WRITE!  1.1120) 

li;0  FORMAT! '  SELECT  !  F 'MOVE, !  P  )3CAN  OR  !T)RACr:  ' ' 

READ! 1. 1110  )ANS 
HLP=3 

IF!ANS.EQ.7:mL  HELP 
IF!ANS.EC!.70)CALL  FMOVE 
IF!ANS.E0.80)CALL  fSCAN 
IFIANS.£Q.84)CALL  TRACK 
IF!ANS.EQ.88)G0T0  190 
GOTO  110 
C 

C  SELECT  TASK  OR  DISPLAY  H.D.E.S.V.X  SCREEN  !riLP=^) 

r 

130  URITE'l.inO) 

1130  FORMAT!'  SELECT  !D)ISPLAY.!E'NTER.!3 'AVE  OR  iV'ELCC:  ' 
PEAD!I.1110.iANS 
HLP-4 

IF!ANS.EQ,72)CALL  HELP 
IF'ANS.E0,48)GOTO  400 
IFfANS.EC(,49)G0T0  600 
IF!ANS.EQ,e3)G0T0  800 
IF<ANS.EO,86)GOTO  900 
iF^4N5.EQ.88)C.OTO  190 
GOTO  130 


c 

C  SELECT  TASK  OR  DISPLAT  H.F.P.T.S.fi.X  SCREEN  (HLP=5) 

140  ORITEl 1,1140) 

1 140  FORHATl '  SELECT  ( F  )HOVE, ( P  )SCAN, ( T IRACL, ( £  )TORE, ( C  lONTINUE:  ' ) 
READ!  1, 1110  )ANS 
HLP=5 

IF(ANS.EQ.72)CALL  HELP 
IF(ANS.EQ.70)CALL  FHOVE 
IF(ANS,EQ.80)CALL  PSCAN 
IF(ANS.EQ.S4)CALL  TRACK 
IF(ANS.EQ.83)G0T0  700 
IF(ANS.EQ.67)GOTO  130 
IF(ANS.EQ.S8)G0T0  190 
GOTO  140 
r 

C  RESTART,  EXIT  OR  CONTINUE 

C 

190  WRITE!  1,1190) 

1190  FORMAT!'  (R)E3TART,  !E)XIT,  OR  (OONTINUE:  ' ' 

READ!  1, 1110  )ANS 
HLP=1 

IF!ANS.EQ.72)CALL  HELP 
IF!ANS.EQ.67)GOTO  130 
IF!ANS.E'].8:)G0T0  100 
IF!ANS.EQ.A9)GOTO  1000 
GOTO  130 
C 

C  ENTER  FILENAME 

C 

300  WRITE!  1,1200) 

1200  FORMAT!  '+ENTER  FILENAME:  ) 

READ!  1,1210)!  FILE!  J).J=1. 11' 

1210  FORMAT!  HAD 

210  WRITE!  1,1220) 

1220  FORMAT!  '+ENTER  !0)LD  OR  !N'EH  FILE:  ‘  * 

READ!  1, 1110  )ANS 
IF!ANS.EQ.79tG0T0  300 
IF(ANS.EQ.78)G0T0  220 
GOTO  210 
C 

C  SETUP  HEADER  FOR  NEW  FILES  OR  CHANCE  OLD  HEADER 


220  CHG=0 

RMAX=0 
FMAX=0 

230  WRITE!  1,  1230) 

1230  FORMAT!  DENTER  name:  ') 

READ!  1,1 23S  )!NAME!J),J'=1,8) 
1235  FORMAT!  8AD 

IF!CHG.EQ.DGOTO  400 
■240  WRITE!  1,1240) 

1240  FORMAT! 'tENTER  DATE:  ') 

READ!  1,1235)!  DATE!  J),  J=l,8) 
IF!CHG.EQ.1)GOTO  400 
250  WRITE!  1,1250) 

1250  FORMAT!  '+ENTER  COMMENTS:  ' ' 
READ!  1. 1260)!  WORDS!  J  >,2=1,46) 
1260  FORMAT!  46A1) 

IF!CHG.EQ.1)G0T0  400 
260  WRITE!  1,1270) 

1270  FORMAT! '+ENTER  FMAX:  ') 

READ!  1, 1275  )FHAX 
1275  FORMAT! 18 ) 

WRITE!  1,1280) 

1280  FORMAT! '  FRAME  TIME'/) 

DO  270  I=1,FMAX 
WRITE!  1,1290  >I 

1300  FORMAT!  F8. 3) 

270  CONTINUE 

GOTO  130 


r 

300 


1310 


READ  OLD  DATA  FROM  DISK 
CALL  0PEN!6,FILE.D' 

READ!  6, 1310  IRMAX.FMA!,;  NAME!  J ).  J=1.E  ),f  DATEi  ',  J=l.  3 
!  WORDS!  J),  2=1,46  I 
FORMAT!  216, 2!  8A1),46A1' 


IF(FfiAX  )370.370,310 
310  DO  320  I=l,FHAX 

READ(  6,1320  )TI«E(I) 

1320  F0RHAT(F8.3) 

320  CONTINUE 

IF(RNAX)370,370, 330 
330  DO  360  I=1.RNAX 

READ(6. 1330)FRN(  I  ),PRT(  I  /.FXl  I  ).Pr(  I  ).PZ(  I  ■. 

1  SAREA(  I  ).SPERIM(  I  ),SEX(  I  ).2EY(  I  ),£«IDEi  I ). STALK  I ), 

2  TAREA(I),TPERIH(I),TLENa) 

1330  F0R«AT<14I8) 

LEN=TLEN( I ) 

IF(LEN  >360,360,340 
340  DO  350  J=1,LEN 

READ(6,1340)TX(I,J),TY(I.J) 

1340  F0RHAT(2I8) 

350  CONTINUE 

360  CONTINUE 

370  ENDFILE  6 

GOTO  130 
C 

C  DISPLAY  FILE  HEADER  AND  CHANGE  AS  REQUIRED 

C 

400  WRITE!  1, 1400  )ESC,  AST,  ( FILE!  J ),  J=l,  1 1 ),!  NANE!  J ),  J=l.3 ). 

1  !  DATE!  J ),  J=1 , 8 ),  RHAX.  FhAX, !  LORDS!  J ),  J=1 . 4i  > 

1400  F0RNAT!1X,2A1/'  FILENAME:  MlAl,'  NAME:  .3A1,  DATE:  ^SAl/ 

1  '  RMAX:  ',18,'  FMAX:  ',18/'  COMMENTS:  ^-oAl 

2  •  FRAME  TIME') 

IF!  FMAX  >420,  420,  405 

405  DO  410  1=1, FMAX 

WRITE!  1,1410)1. TIME!  I) 

1410  F0RMAT!I4,4X.F3.3/ 

410  CONTINUE 

420  WRITE!  1,1420) 

1420  FORMAT!'  CHANGE  !N  )A.i£.!  D 'ATE.IC  lOMilENTS  OR  'F'SAME  TIME’  ') 

READ!  1. 1110  )ANS 
CHG=1 
HLP=7 

IF!ANS.EQ.72>CALL  HELP 
IF!ANS.EQ.78'G0T0  220 
IF!ANS.£Q.68>GOTO  240 
IF!ANS.EQ.67)G0T0  250 
IF!ANS.£a.70)G0T0  422 
GOTO  430 

422  WRITE!  1,  1424  ) 

1424  FORMAT! '4ENTER  FRAME  NO:  ') 

READ!  1,1275  )I 
IF! I -FMAX >428,428, 426 
426  FMAX=FMAX4l 

TIME!  I  )=0 

423  WRITE!  1, 1430  )TIME!  I) 

■1430  FORMAT! '40LD  FRAME  TIME:  ',F8.3/'  ENTER  NEW  FRAME  TIME:  ') 

READ!  1,  1300  )TIME!I) 

GOTO  400 
C 

C  ENTER  FRAME  AND  PARTICLE  FOR  DISPLAY 

C 

430  WRITE!  1,1440) 

1440  FORMAT!  '+ENTER  FRAME  NUMBER  OR  !CR)  FOR  ALL  FfAMES:  ) 

READ!  1, 1275  )IFPH 
WRITE!  1,  1450) 

1450  FORMAT!  '+ENTER  PARTICLE  NUMBER  OR  !CP)  FOR  ALL  PARTICLES:  ' > 
READ!  1,  1275  )IPRT 
IF!  IFRN  >432,440,432 
432  IF!  IPRT  >540,510, 540 

440  IF!  IPRT >480,450, 480 

r 

•;  DISPLAY  ALL  BATA  IN  FILE  BY  FRAME  AND  PARTICLE 

r 

450  WRITE!  1,1460  >ESC, AST, ! FILE! J  >, J=l,  11 > 

1-60  FORMAT!  IX.  2A1/'  HLENAME:  ',11A1/ 

1  '  — FRAME-PAPTCLE-x-AXIS-Y  AXIS-Z  AXIS', 

2  '-SAREA’-SPERIM—SWIDE— STALL' ' 

00  470  1=1, RMAX 

460  WRITE!  1.  1470  'FRM!  I  ',RRT(  1 1  p<(  I  \F-ff  I  ,>,P:!  I ', 

1  SAf EA'  I  >. SPERIM'  I  >.  3WIDE'  I  V'S'^ALL!  I ' 
l-’O  FORMAT! ’13* 

♦■’0  CONTINUE 

GOTO  130 


DISPLAY  PARTICLE  DATA  BY  FRAME 

WRITE!  1,1480  )ESC,  AST,  ( FILE!  J ).  J=l,  11 ),  IPRT 
FORMAT!  IX, 2A1/'  hLEMAHE:  ',11A1,"  PARTICLE:  ',13/ 

'  — FRAME-X  AXIS-Y  AXIS--Z  AXIS', 

'  - — 5CMX- — SCNY — SAREA--SPERIM—  EWIDE— STALL ' ' 

DO  500  I=1,RMAX 
IF!PRTII)-IPRT)500,490,500 
WRITE!  1,1490  )FRM!  I  \  PX!  I ),  PY!  I ),  PZ!  I ), 

SEX!  I ), SEY!  I ), SAREA!  I ),  SPERIM!  I  ),£WIDE!  I ),  STALL!  I ) 
FORMAT!  1018) 

CONTINUE 
GOTO  130 

DISPLAY  FRAME  DATA  BY  PARTICLE 

WRITE!  1 , 1500  )ESC,  AST, !  FILE!  J ),  J=l.  11 ),  IFRM 
FORMAT! IX, 2A1/'  FILENAME:  '.llAl/'  FPAnE:  ',18/ 

'  -PARTCLE-X  AXIS-Y  AXI3-Z  AXIS', 

'— SAREA--SPERIM — 5WIDE— STALL— TAREA-TPERIM') 

DO  530  I=1,RMAX 
IF!  FRM!I)-IFRM)530,520.530 
WRITE!  1, 1490  )PRT!  I  ),PX!  I  ),PY!  I  ),P2!  I ), 

SAREA!  I ),  SPERIM!  I ),  SWIDE!  I ),  STALL!  I ), 
TAREA!I),TPERIM!I) 

CONTINUE 
GOTO  130 

DISPLAY  EDGE  DATA  FOR  SELECTED  PARTICLE  AND  FRAME 

DO  560  I=1,RMAX 
IF!FRN!n-IFRM  )560.550,5i0 
IF!PRT!n-IPRT  >560,570,540 
CONTINUE 
GOTO  130 

WRITE!  1, 1510  )ESC,AST,  (FILE!  J !,  J=l,  11 ), 

IFRM,  IPRT,  I,  TAREA!  I ),  TFERIM!  I  ),TLEN(  l' 

FORMAT!  IX.2A1/ 

'  FILENAME:  '.llAl/'  FRAME:  M8,"  PARTICLE:  ',18/ 

'  RECORD:  ',18/'  TAREA:  MS,"  TFERIM:  ',IS"  TLEN:  ',IS 
'  -POINT-’X  AXIS-Y  AXIS') 

LEN=TLEN!  I ) 

IF! LEH >130,130,530 
DO  590  J=1,LEN 

WRITE!  1,1520)J,TX!I,J),TY(I,J) 

FORMAT! 318) 

CONTINUE 
GOTO  130 

ENTER  FRAME  AND  PARTICLE  IDENTIFICATION 
WRITE!  1,1600) 

FORMAT! '+ENTER  FRAME  NO:  ') 

READ! 1, 1275 )FRAME 
IF!  FRAME  )600, 600, 610 
IF!  FRAME-9 >620, 620, 600 
WRITE!  1,1610) 

FORMAT! 'tENTER  PARTICLE  NO:  '' 

READ!  1,1275  )PART 
IF!  PART >620,620, 630 
IF!  PART-9  >640,640, 420 
DO  650  I=1,RMAX 
IF!FRM!I)-FRAME)630,645,i50 
IF!  PRT!  I  )-PART  >650,660, 650 
CONTINUE 

RESET  FRAME  AND  PARTICLE  DATA  TO  ZERO 


TH( J )=0  ' 

TH(J)=0 
CONTINUE 
GOTO  140 

RESET  TO  RECORD  DATA 

FX=PX<  I ) 

FY=PY( I ) 

FZ=PZ(  I ) 

SX=SEX(  I ) 

SY=SEi(I) 

SA=SAREA( I ) 

SP=SPERIM( I ) 

SU=SUIDE( I ) 

SH=STALL( I ) 

TA=TAREA(I) 

TP=TPERIN( I ) 

TL=TLEN( I ) 

IF(  TL)465.670.46C 
DO  ft70  J=1,TL 
TU(J)=TX(I,J) 

TH(J)=TY(I,J) 

CONTINUE 
GOTO  140 

SET  RECORD  DATA  FOR  SELECTED  FPAHE  AIID  PARTICLE 

IF(RHAX  )730.730.710 
DO  730  I=1.RHAX 
IF(FRM(I)-FPAME  >730.720,730 
IF(  PRT(I)-PART)730,  740.  730 
CONTINUE 
RMAX=RHAX+1 
I=RHAX 

FRN(  I  )=FRAhE 
PRT(  I  )=PART 
PX(  I  '=FX 
PY(  I  >=FY 
PZ(  I  )=FZ 
SEX(  I  )=SX 
SEY(I)=SY 
SARf A(  I  )=SA 
SFERl«(n=SF 
SUIDE(  I  )=SI4 
STALL!  I  '=SH 
TAREA(  I  )=TA 
TPERII1(I)=TP 
TLEH(  I  )-TL 
IF(TL  )750,7i0,750 
DO  760  J=1.TL 
TX(I,J)=TU(  J) 

TY(i;j)=TH(  J) 

CONTINUE 
GOTO  130 

SAVE  DATA  ON  DISK 
CALL  0PEN(6,FILE.l) 

URITE! 6, 1810  )RMAX, FMAX, (  NAME!  J  >,  J=I , 8 ),( DATE!  J J^l. S  i. 
!W0RDS!J),J=1,46) 

F0RMAT!1X,:I8,2!3A1),46A1) 

IF! FNAX >860.860,810 
DO  820  1=1,  FMAX 
URITE!6,i820)TIME!I) 

FORMAT!  1X,F6. 3) 

CONTINUE 

IF!RMAX)860,860,830 
DO  860  I=1,RNAX 

WRITE!  6,  1830  >FRM!  I  ).FRT!  I  ).PX!  I  ;,PY!  I  ).f  Z!  I ', 

SAREA!  I  ),SPERIM!  I  ).SEX!  I  ).SEY!  I  .'.SMIDE!  I ). STALL!  I  >. 
TAREA!I),TPERIMII),TLEN!I) 

FORMAT!  IX, 1418) 

LEN=TLEN! I ) 

IF!LEN  >860,860,840 
DO  850  J=1,LEN 
WRITE!  6. 1840  )TX!I,J),TY!I.J) 

FORMAT! IX, 218) 

CONTINUE 


ENDFILE  i 
GOTO  110 

COMPUTE  DRIFT  VELOCITY  FOR  SELECTED  PARTICLES  AND  FRAMES 

WRITEd,  1900  )ESC,  AST 
F0RMAT(1X.2A1, 

■  DRIFT  VELOCITY  AND  AVERAGE  PARTICLE  SIZE  COMPUTATION'/ 
'  ENTER  FIRST  FRAME:  ' ) 

READ(1.1910)FR1 
FORMAT!  18) 

«RITE!1.19:0' 

FORMAT! '+ENTER  SECOND  FRAME:  ') 

READ!  1. 1910  )FR: 

DT=TIME!FR2)-TIHE!FR1) 

WRITE!  1. 1930  )DT 

FORMAT! '+TIME  DIFFERENCE:  '.'3.2,  / 

'  PARTICLE— X  VEL — Y  VEL— Z  VEL'. 

— AREA — FERIM — WIDTH-HEIGHT ' ' 

IFIDT)902.130.902 
IF!  RMAX  )130,  130. 905 
DO  950  1=1. RMAX 
IF!FRM!I)-FR1  )950.?10.?30 
DO  940  J=1.F:MAX 
IF!  FRM!J)-FR2  )940.?:0.?40 
IF!PRT(I)-PRT!J))9-t0.dO.’:40 
DX=PX(  J  '-PX!  I ) 

DY=PY!J)-PY(I) 

DZ=PZ(J)-PZ!n 

XVEL=DX/DT 

YV£L=DY/DT 

ZVEL=DZ/DT 

AREA=(  SAREA!  I  )+SAREA<  J )  :-,'2 
PERIM=(SPEPIM<I)^ZF'ERIM!J))/: 
HIDTH=!SgiD£li:'+5Un<[(jn/2 
HEIGHT*!  STALL!  I  INSTALL!  J}'/; 

WRITE!  1.  1940  )PRT!  1 /';EL.  YVZL.ZVEL, 

area.perim.width.height 

FORMAT!  IX,  18. 7F8. 3) 

CONTINUE 
CONTINUE 
GOTO  130 

EXIT  FROM  PROGRAM  AND  RETURN  TO  SYSTEM 
CONTINUE 


c 

c 

c 

c 

c 

c 


c 

1000 

c 

c 

c 

100 

1010 

1 

T 

3 

4 

j 

6 

110 

1020 

1 

■) 

3 

4 

5 

6 

7 

8 
9 
A 

C 

C 

C 

200 

1200 

1 

1 

3 

4 

5 

6 

7 

8 
9 

C 

C 

C 

300 

1300 

1 


4 

C 

C 

r 

400 

1400 

1 

3 

4 

C 

I 

7 


SUBPROGRAM  TO  HELP  OPERATOR  USE  THIS  SYSTEM  BY  DISPLAYING 
ADDITIONAL  INFORMATION  AND  INSTRUCTIONS  ON  THE  CRT  UHEN 

call;d  from  different  parts  of  the  program 

REVISION  0.12,  28  AUC  1982 

SUBROUTINE  HELP 

LOGICAL  ESC. AST 
INTEGER  HLP 
COMMON  ESC, AST, HLP 

HRITEI I  1000  )ESC  AST 

FORMAT!  ix,2Al,  ADDITIONAL  INFORMATION  AND  INSTRUCTIONS'-') 
GOTO!  100, 200, 300, 400, 500, 600, 700,200 ),HLP 

OVERVIEW  SCREEN 

WRITE!  1.1010 'ESC,  AST 
FORMAT!  IX,  2A1/ 

'  - HOLOGRAPHIC  IMAGE  ANALYSIS  SYSTEM . . 

'  Written  by'/ 

AZTEC  COMPUTER  ENGINEERING'/ 

'  St.  Pitersbur-i,  Floiida'.' 

'  28  AUGUST  I?92'/ 


WRITE!  1,1020) 

FORMAT!'  THE  PROGRAM  WILL  PERFORM  THE  FOLLCWING  FlvCTIONS: 

'  1.  IH)ELP  DISPLAYS  MORE  DETAIL  AT  ANY  STEP  IN  THE  PROGRAM  / 

'  2.  (L)OADS  EXISTING  DATA  FROM  DISK  TO  IIEnORY'/ 

■■  3.  !D)I3PLAY  PRINTS  DATA  BASE  BY  FRAME  OR  PARTICLE  ON  CRT', 

'  4.  !E  INTER  ADDS  NEW  DATA  OR:  CHANGES  CLD  DATA  IN  pEMORY' 

'  5,  1S)AVE  WRITES  COMPLETE  DATA  BA:SE  FROM  MEMORY  TO  DISK'  ' 

'  6.  (F)RAME  MOVES  STAGE  TO  REFERENCE  OR  FARTICLE  POSITION'/ 

'  7.  !P)ARTICLE  SCANS  IMAGE  TO  FIND  PARTICLE  SIZE 
'  7.  !T)RACK  TRACES  IMAGE  TO  FIND  PARTICLE  SIZE  AND  EDGE 
'  8.  IVIELOCITY  COMPUTES  DRIFT  VALUES  FOR  SELECTED  f ARTICLES 
'  10.  E(X)II  LEAVES  PROGRAM  AND  RETURNS  TO  OPERATING  SYSTEM  ,> 
RETURN 

LOAD  AND  TEST  SCREEN 

WRITE!  1, 1200  )£SC, AST 
FORMAT!  IX, 2A1/ 

'  THE  FOLLOWING  FUNCTIONS  MAY  BE  INVOKED  AT  THIS  LEVEL:'/-' 

'  1.  (H)£LP  DISPLAYS  THIS  SCREEN'// 

'  2.  IDOAD  READS  NEW  DATA  FROM  DISK  OVER  OLD  DATA  IN  MEMORY'// 
'  3.  !T)EST  ENABLES  THE  STAGE,  TUSCAN  AND  TVTRACK  FUNCTIONS'/ 

'  TO  BE  TESTED  BEFORE  ANY  DISK  FILES  ARE  LOADED  IN  MEMORY'// 
'  4.  E!X)IT  RETURNS  OPERATION  TO  CDOS','/ 

'  REMEMBER!"'// 

'  SAVE  OLD  DATA  BEFORE  LOADING  NEW  DATA  CR'/ 

EXITING  TO  OPERATING  SYSTEM!!  //) 

RETURN 

FRAME,  SCAN  AND  TRACK  SCREEN 

WRITE!  1,1300  )ESC,AST 
FORMAT!  IX, 2A1/ 

'  THE  FOLLOWING  FUNCTIONS  MAY  BE  TESTED  AT  THIS  LEVEL:'/-' 

'  1.  (F)RAME  ENABLES  XYZ  STAGE  TO  BE  TESTED  WITH  JOYSTICK'.'/ 

'  2.  !P PARTICLE  ENABLES  IMAGE  TO  BE  SCANNED  WITH  TV  CURSOR'-'/ 

'  3.  (T)RACK  ENABLES  PARTICLE  EDGE  TO  li  TRACED  BY  PROGRAM'//  ) 
RETURN 

DISPLAY,  ENTER,  SAVE  AND  VELOCITY  SCFEEN 


WRITE! 1, 1400  )ESC, AST 
FORMAT!  lx, 2A1/ 

'  THE  FOLLOWING  FUNCTIONS  MAY  BE  INVOKED  AT  THIS  LEVEL:' 

'  1.  (D)ISPLAY  PRINTS  POSITIONAL  DATA  li  FRAME  OR  PAPIICLE  ' 

'  AND  EDGE  DATA  UHEN  BOTH  FRAME  AND  PAPTICLE  ARE  ENTERED  // 
'  2.  !E  INTER  ENABLES  FRAME  AND  PARTICLE  DATA  FROM  XYZ  STAGE, 

'  TV  SCAN  AND  TRACK  FUNCTIONS  TO  EE  IIENTIFIED  AND  STORED  - 
'  3.  !S)AV£  WRITES  ALL  DATA  II!  MEMORY  TO  DISK  LNJER  FILENAME', 
'  4.  (V)ELOCITY  COMPUTES  DRIFT  VALUES  FOR  SELECTED  PARTICLES  , 
RETURN 


C 

c 

r 


FRAME.  SCAN,  TRACK,  SET  AND  RESET  SCREEN 


TO  MOVE  STAGE  IN  XTZ  AXIS 
REVISION  0.28,  2  APR  1983 

?•  <R)RA«E  REFERENCE  (MOVE  X.Y.Z  TO  0.0.0) 

8.  nOVE  STAGE  TO  (P)ARTICLE  POSITION  ChOVE  X.Y.Z  TO  FX.PY.FZ) 
ENTER  COMMANDS  FROM  (  J)OYSTICK 

IN  X  AND  Y  AXIS 

2.  MOVE  STAGE  IN  X  AND  Z  AXIS 

3.  SET  PARTICLE  POSITION  AND  RETURN 

4.  RETURN 

D.  ENTER  COMMANDS  FROM  (K)EY6OAR0 

1.  INPUT  (X).(Y).(2)  AXIS 

2.  (R)ESET  FRAME  REFERENCE  f RESET  X.Y.Z  TO  0.0.0) 

3.  (S)ET  PARTICLE  POSITION  <3ET  FX.FY.FZ  TO  X.Y.Z) 

4.  RETURN 

E.  RETURN 

SUBROUTINE  FMOVE 
LOGICAL  E3C.AST.ANS 

INTEGER  XYZM. XOUT. TOUT. ZQUT. K JS,  JSK.  JSX. JSY.  JSZ 

INTEGER  HLP, FX. FY. FZ. SX. 3Y. SA, SP. S«. SH. TA. TP.  TL,TU(  100 ). TW  1 00 ) 

EXTERNAL  J3INP.STAGE 

COMMON  ESC. AST. HLP, FX, FY. FZ, SX, 3Y, SA. SP. SU. SH. TA. TP. TL.  TU.  TH 

ENTRY  POINT 

XYZM=0 

CALL  TVSCN(XYZH) 

X0UT=0 

Y0UT=0 

Z0UT=0 

JSK=0 

JSX=0 

JSY=0 

KJS=0 

SELECT  MODE  FOR  MOVING  STAGE 
URITE(  1.1100) 

FORMAT! '+MOVE  TO  (F)RAHE  REFERENCE. (DARTICLE  POSITION.'. 
'(J)OYSTICK  OR  (K)£YBOARD  ENTRY:  -  ) 

READ!  1, 1110  )ANS 
FORMAT! lAl ) 

HLR=11 

IF!ANS.EQ.72)CALL  HELP 
IF!ANS.EQ.70)GOTO  200 
IF!ANS.EQ.80)G0T0  300 
IF!ANS.EQ.74)G0T0  400 
IF!ANS.EQ.75)G0T0  500 
GOTO  600 

MOVE  STAGE  TO  FRAME  REFERENCE 

XYZM=1 

X0UT=0 

youT=o 

Z0UT=0 

CALL  STAGE!  XYZM. XOUT. YOUT) 

WRITE!  1,1200) 

FORMAT! ^MOVE  STAGE  TO  FRAME  REFERENCE') 

GOTO  100 

MOVE  STAGE  TO  PARTICLE  POSITION 

XYZN=1 

XOUT=FX 

YOUT=FY 

ZOUT=FZ 

CALL  STAGE! XYZM. XOUT. YOUT) 

WRITE!  1,1300) 

FORMAT' '♦MOVE  STAGE  TO  PARTICLE  POSITION') 

GOTO  100 

ENTER  COMMANDS  WITH  JOYSTICK 

WRITE!  1. 1400  )E3C, AST 
FORMAT!  IX, 2AI, 

'  . . JOYSTICK  XYZ  MODES . . ' 

'  1.  MOVE  STAGE  RIGHT  !X)  OR  UP  lY)', 


2.  HOVE  STAGE  RIGHT  (X)  OR  FORHARD  (Z)'/ 

'  3.  SET  POSITION  '< 

'  4.  CONTINUE'/ 

jg-“-NODE - JSX - JSY . JSZ - XOUT- — TOUT - ZOUT ' 

JSX=6 

JSY=0 

JSZ=0 

KJS=0 

«RITE(  1,1410  )K  JS.  JSX,  JSY.  JSZ.  XOUT.  YOUT.ZCUT 
F0RHAT('+',7I8) 

CALL  JSIHP(  JSK.JSX.JSY) 

IF(  JSK >430,430, 420 
IF(JSK-9  >440.600.600 
IF(KJS>410.410,450 
KJS^JSK 

GOTO(  460. 470, 470, 560, 560, 560, 560, 100  >. K JS 

X  AND  Y  AXIS 

XYZH=1 

XOUT=X0UT+JSX/10 
Y0UT=Y0UT+JSY/10 
CALL  STAGE! XYZN. XOUT. TOUT) 

GOTO  410 

X  AND  Z  AXIS 

XYZH=1 

JSZ=JSY 

XOUT=XOUT+JSX/10 
ZOUT=Z0UT+JS2/10 
CALL  STAGE!  XYZN. XOUT. TOUT > 

GOTO  410 

ENTER  CCHNANDS  FROM  KEYBOARD 
XYZM=1 

WRITE!  1, 1500 > 

FORMAT! '♦ENTER  !X>.!  ,;,!Z)  VALUES. !R>ESET  REFERENCE  OR'. 
'!S)ET  position:  '  > 

READ!1.1110>ANS 

HLP=12 

IF!ANS.EQ.72>CALL  HELP 
IF!ANS.£Q.88>GOTO  510 
IF(ANS.EQ.8?>G0T0  520 
IF!ANS.EQ.90>G0T0  530 
IF!ANS.EQ,82>G0T0  550 
IF!ANS.EQ.83>GOTO  560 
GOTO  100 

WRITE! 1, 1510 > 

FORMAT!  '+X  VALUE:  ' > 

READ!  1. 1520  >XOUT 
FORMAT!  18 > 

GOTO  540 
WRITE!  1,1530) 

FORMAT! '♦Y  VALUE:  '> 

READ!  1, 1520  )YOUT 
GOTO  540 
WRITE!  1,1540) 

FORMAT! '>3  VALUE:  ') 

READ!  1,  1520  >ZOUT 

CALL  STAGE! XYZM, XOUT. TOUT) 

GOTO  500 

XYZM=0 

CALL  STAGE! XYZM) 

GOTO  500 

FX=XOUT 
FY=YOUT 
FZ=ZOUT 
GOTO  lOO 


CONTINUE 

RETURN 

END 


STAGE  SUBROUTINE 
REVISION  0.20,  27  JULY  1982 

THIS  ROUTINE  HOVES  THE  XY2  STAGE  TO  THE  DESIGNATED 
POSITION  AND  RETURNS. 

STAGE!  XYZNCHL2, XOUTEDEJ. YOUTEBCJ  ),Z0UTCBC+:3 


ENTRY 

STAGE:  LD 
ADD 
JP 
LD 
OUT 
RET 

^ART:  PUSH 
POP 
INC 
INC 
PUSH 
PUSH 
PUSH 
LD 
LD 
LD 
LD 

AXIS:  LD 


STAGE 

A,(HL) 

A,0 

NZ, START 
A.O 

(§4H>,A 


BC 

HL 

HL 

HL 

HL 

BC 

DE 

A,  208 
(ACODE),A 
A, 3 

(ACOUNT),A 

A,  5 

(tiCOUNT),A 

C, 192 

hL 

E,(HL) 

HL 

D, (HL) 

dLhl 

IX, DECADE 

E, (IX) 

IX 

D, (IX) 

IX 

HL,DE 

P,NEG 

HL,D£ 

B. O 
POS 

6,3 

DE,HL 

HL,DE 

A 

A 

E, (IX) 

IX 

D,(IX) 

IX 

A 

HL,DE 

P,LOOP 

A 

HL,DE 

A 

A,B 

A,C 

(84H),A 

A,(DCOUNT) 

Z, OUTPUT 


;GET  FLAG 
iCHECK  ZERO  FLAG 
: START  XFER  IF  FLAG  OK 
; RESET  A  PEG 
:AND  OUTPUT 
; RETURN 

;GET  Y  AXIS  ADDRESS 
:AND  PUT  IT  IN  HL  REG 
SINCRENENT  ONCE 
(TWICE 

(SAVE  Z  AXIS  ADDRESS 
(SAVE  Y  AXIS  ADDRESS 
(SAVE  X  AXIS  ADDRESS 
( INIT  AXIS  CODE 
(AND  SAVE 
(INIT  AXIS  COUNT 
(AND  SAVE 
(INIT  DIGIT  COUNT 
(AND  SAVE 
(INIT  DIGIT  CODE 
RECALL  AXIS  ADDRESS 
(GET  AXIS  VALUE  (LOW  BYTE) 

((HIGH  BYTE) 

(MOVE  VALUE  TO  HL 

(INIT  DECADE  ADDRESS 

(LOAD  HIGHEST  DECADE  (LOU  BYTE) 

((HIGH  BYTE) 

(SUBTRACT  HIGHEST  DECADE  FRO!i  VALUE 
(JUrtP  IF  VALUE  IS  NEGATIVE 
ADO  DACK  HIGHEST  DECADE 
(LOAD  POSITIVE  FLAG 
(VALUE  IS  positive 

(LOAD  NEGATIVE  FLAG  INTO  B  REG 
(EXCHANGE  REGS 
(CONPLEHENT  VALUE 
(RESET  OUTPUT  WORD  TO  ZERO 
(CLEAR  CARRY 

(LOAD  NEXT  DECADE  (LOU  BYTE) 

((HIGH  BYTE) 

(INCREHENT  DIGIT  VALUE 
(SUBTRACT  NEXT  DECADE  FROM  VALUE 
(DO  AGAIN  IF  STILL  FG3IT1VE 
(CLEAR  CARRY 
(TOO  HUGH,  ADD  IT  BACK 

(ADD  SIGN  TO  OUTPUT  WORD 
(ADD  DIGIT  CODE  TO  OUTPUT  UORD 
(OUTPUT 

(GET  DIGIT  COUNT 

(AND  DECREMENT 

(OUTPUT  AXIS  CODE  IF  DONE 


LD 

(OCOUNT).A 

(SAVE  DIGIT  COUNT 

t! 

w 

(RESET  SIGN  FLAG  TO 

f 

IK' 

LD 

A,(ACODE) 

(LOAD  AXIS  CODE 

OUT 

(84H),A 

(OUTPUT 

ADD 

A.IOH 

( INC  AXIS  CODE 

LD 

(ACODE),A 

(AND  SAVE 

kcADE:  DM  8000H,  10000, 1000, 100, 10,1 

flCODE:  DS  2 

ACOUHT:  DS  2 

DCOUNT:  D3  2 


TYPE  JSIMP.ZSO 
:  JOYSTICK  SUBROUTINE 

;  VERSION  0.21,  le  JULY  1962 

;  THIS  ROUTINE  INPUTS  UNIPOLAR  DATA  FROM  THE  MODE  SWITCH.  X  AX 

AND  Y  AXIS  THEN  RETURNS  BIPOLAR  VALUES  TO  CALLING  PROGRAM 


JSINP: 


JSXY! 


J30UT; 


RETURN! 


ENTRY 

PUSH 

PUSH 

PUSH 

IN 

LD 

LD 

SUE 

POP 

LD 

LD 

LD 

IN 

LD 

LD 

LD 

SEC 

IP 

EX 

LD 

SEC 

EX 

POP 

LD 

INC 

LD 

DEC 

JP 

INC 

JP 

RET 

END 


JSINP 

EC 

DE 

HL 

A. ( 18H ) 
E.  A 
A.  255 
E 

HL 

( HL) . A 
E.  2 

C.  I'-H 
A,  (C  ) 

E.  A 

D.  0 

HL. 12B 
HL.  DE 
P. JSOUT 
DE.  HL 
DE. 255 
HL.  DE 
DE.  HL 
HL 

( HL) . E 
HL 

(HL) . D 
B 

Z. RETURN 
C 

JSXY 


SAVE  JSV  ADDRESS 
SAVE  JSX  ADDRESS 
SAVE  J5K  ADDRESS 
INPUT  JS  SWITCH 
MOVE  IT  TO  B 
LOAD  A  WITH  255 
SUE  JS  SWITCH 
GET  JSK  ADDRESS 
STORE  MODE 
INIT  AXIS  COUNT 
INIT  INPUT  PORT 
INPUT  JS  VALUE 
LOAD  E  WITH  JS  VALUE 

LOAD  HL  WITH  12B 
SUBTRACT  JS  VALUE  FROM  123 
JUMP  IF  POSITIVE  VALUE 
FIND  NEGATIVE  VALUE 


;GET  X  OR  Y  AXIS  ADDRESS 
.-LOAD  BIPOLAR  VALUE  FOE  RETURN 


;MOVE  TO  NEXT  AXIS 
; RETURN  IF  DONE 
;MUST  EE  Y  AXIS 


PSCAN  SUBROUTINE  TO  INPUT  PARTICLE  SCAN  DATA 
REVISION  0.28,  7  APR  1983 

SUBROUTINE  PSCAN 


r 


c 

c 

c 


c 


c 

c 

c 


c 

c 

c 

r 


1000 


c 

c 

c 

10 

20 

30 

40 

50 

C 

C 

C 

100 

no 

120 


c 

c 

c 

200 

210 

220 


LOGICAL  ESC.  AST, ANSWER 

INTEGER  KJS, JSK, JSX , JSY. THRESH. VIDEO 

INTEGER  XOUT, TOUT, LOUT, TOUT, VINP, AINP,  PINP 

INTEGER  HLP.PX.FY,F2,SX,SY,SA,£P.SW,£H,TA.TP,TL.TU(  100).TH(  100) 

COMMON  ESC, AST, HLP, FX, FY, F2, EX, 3Y, Sft, SP. iW, SH, TA, TP, TL, TU. TH 

ZERO  CURSOR  OH  TV  MONITOR 


JSK=0 

JSX=0 

JSY=0 

KJS=0 

LR=0 

IPRINT=0 

MR=10 

NR=10 

SX=123 

SY=128 

SW=1 

SH=1 

X0UT=i:8 

YOUT=128 

L0UT=130 

TOUTsO 

VINP=0 

AINP=0 

PINP=0 

VIDEO=0 

THRESH=0 

CALL  TVSCNITOUTI 

DISPLAY  JOYSTICi:  SELECTIONS 


WRITE!  1,1000  'ESC, AST 
FORMAT! IX, 2A1, 

1  ' - JOYSTICK  PARTICLE  SCAN  MODES 


n;- 


?  '  r.  5?  PARTICLE  - 

4  '  3.  MOVE  CENTER  OF  CURSOR  SCAN'/ 

5  '  4.  ENTER  COMMANDS  WITH  KEYBOARD'/ 

A  '  -JSMODE--X  CNTR--Y  CNTR—WIDTH--I-EIGHT— VIDEO'.' 


1 


FIND  PARTICLE  CENTER,  SIZE  AND  SHAPE 


KJS=5 

CALL  JSINP!  JSK,JSX,JSY) 

IF!  JSK)20,20,30 

IF!KJS)10,10.50 

IF!JSK-9  >40,900,900 

KJS=JSK 

LR=0 

GOTO!  100,200, 200, 300,  .^00, 300, 300, 400 ',  K  JS 
MOVE  CURSOR  WITH  JOYSTICK 


IF!LR)110, 110,120 

XOUT=SX 

YOUT=SY 

LRM 

XOUT=XOUT+JSX/MR 

YOUT=YOUT-JSY/Mfi 

L0UT=Y0UT+1 

T0UT=1 

CALL  TVSCN! XOUT, 'JIMP ) 

VIDEO=VINP 

GOTO  700 


FIND  PARTICLE  SCAN  SIZE 


SH=SH+JSX/HR 
IF!SH  '210,210,220 
£W=1 

SH=SH+JSY/HR 
IF!  SH  >230,230, 240 


220  SH=1 

240  X0UT=SX-SU 

yOLiT=SY-SH 
L0iJT=St+3H 
T0LIT=1 

CALL  TVSCN(XOUT,VINP) 

X0UT=SX 

CALL  TVSCN(XObi'.VIHP) 
XOUT=SX+SW 

CALL  TVSCN(XOUT,VINP) 

GOTO  700 
C 

C  FIND  PARTICLE  SCAN  CENTER 

C 

300  SX=SX+JSX/NR 

SY=SY-JSY/NR 
YOUT=SY-SH 
LOUT=SY+SH 
T0'JT=1 
XOUT=SX-S« 

CALL  TOSCNvXOUT.VINP) 

X0UT=SX 

CALL  TVSCN(XQUT,yiNP) 
X0UT=SX+S« 

CALL  TVSCN(XOUT,giNP) 

GOTO  700 

c 

C  ENTER  COMMANDS  WITH  KEYBOARD 

C 

400  «RITE(  1.1400) 

1400  FORMAT!  IX, 


1  '  - KEYBOARD  COMMANDS - '/ 

2  '  1.  ENTER  THRESHOLD  AND  SCAN  IMAGE'/ 

3  '  2.  REPEAT  SCAN  WITH  SAME  THRESHOLD'-' 

4  '  3.  JOYSTICK  CONTROL  WITH  IISPLAr  -' 

5  '  4.  JOYSTICK  CONTROL  WITHOUT  DISPLAY'/ 

6  '  5.  SET  VALUES  AND  RETURN'-') 


WRITE!  1.1410) 

1410  FORMAT!  IX. 'ENTER  SELECTION!  '' 

READ!  1. 1420  )ISEL 
1420  FORMAT! 14) 

KJS=0 

JKS=0 

IFIISEL  )400.4C0,410 
410  IF!ISEL-4)420.  420.400 

420  GOTO!  500, 550, 600, 650, ?00 ), I3EL 

C 

C  ENTER  THRESHOLD 

C 

500  WRITE!  1,1500) 

1500  FORMAT! '  ENTER  THRESHOLD:  ' ' 

READ!!.  1420  )THRF5H 
IF!  THRESH  )500,500,510 
510  IF! THRESH-255 >800.800. 500 

550  IF!  THRESH  >500,500, 800 

C  SET/RESET  PRINT  FLAG 

C 

600  IPRINT=0 

GOTO  5 
C 

650  IPRINT=1 

GOTO  5 

r 

C  PRINT  PARAMETERS  IF  FLAG  IS  SET 

C 

700  IF!IPRINT)710.710,10 

710  WRITE!  l.l 710  )KJS,SX,SY,5y.SH, VIDEO 

1710  FORMAT! '  +  ',618) 

TOUT^O 

CALL  TVSCH!XOUT) 

GOTO  10 
r 

750  WRITE!  1,1750  ) 

1750  FORMAT! ix. 

I  '  -KB  MODE— X  CTR— Y  CTR— -AREA— PERIM-THRESH 


1760 

C 

C 

C 

800 


810 

820 

330 

C 

C 

C 

?00 


“RITE!  M760  )ISEL,SX,SY, £A,  Sr,  THRESH 
FWHAT(6I8) 

GOTO  600 

SCAN  PARTICLE  FOR  SIZE  AND  SHAPE 

XOUT=SX-SH 

YOUT=SY-SH 

LQUT=SY»SH 

XLII1=SX+3y 

T0UT=THRESH 

SA=0 

SP=0 

LASTA=0 

CALL  TVSCN(XOUT,W) 
SA=SA+(AINP»0.79) 

SP=3P+PINP+(  IABS(  LASTA-AINP )t0.79 ) 

LASTA=AINP 

X0UT=X0UT+1 

IF(  XOUT-XLIH)810,830,330 

KJS=0 

GOTO  750 

RETURN  TO  CALLING  FPOGRAM 

CONTINUE 

RETURN 

END 


A. 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

10 

11 


12 


20 

1000 

1 

2 

3 

4 

5 
i 
7 


30 

1010 

40 

SO 

60 

70 

80 

C 

C 

C 

100 


c 

c 

c 

200 


TRACK. FOR  SUBPROGRAH 
REVISION  0.34,  7  FEB  83 

THIS  PROGRAH  PROVIDES  FOR  AUTO  TRACKING  A  VIDEO  IHACE. 

SAVING  THE  EDGE  COORDINATES  OF  THE  VIDEO  IMAGE, 

COMPUTING  THE  PERIMETER  AND  AREA  OF  THE  IMAGE  AND 
DISPLAYING  THE  PROCESSED  PERIMETER  ON  TOP  OF  THE  IMAGE. 

SUBROUTINE  TRACK 

LOGICAL  ESC, AST 
INTEGER  JSK,JSX,JSY 

INTEGER  XOUT,YOUT,LOUT.TOUT,VIHP.AINP,PINP 
INTEGER  SCAM^  SCART,  START,  AW,  DU,  DH,  DX(  64 ),  Dr ( 64  > 

INTEGER  NSCAH, TOLER, POINT, XSCAHl 32\ YSCAH( 32 ) 

INTEGER  AREA, PERIM,PTSI!E, VIDEO 

INTEGER  HLP, FX,FY, FZiSX, SY, SA, SP, SU, SH.TA JP,TL. TWl  100 TH(  100 ) 
COMMON  esc,ast,hlp,N,fy,fz,sx,sy,sa,sf,Xsh,ta,tf,tl,tu,th 

DISPLAY  JOYSTICK  SELECTIONS 


WRITE!  1,11) 

FQRHAT! '  ENTER  FY/FX:  ' ) 

READ<1.12)FH 
FORMAT!  F8. 3) 

KJS=1 

JSK=0 

JSX=0 

JSY=0 

X0UT=128 

Y0UT=123 

VIDE0=0 

TL=0 

PTSIZE=TL 

PER1M=0 

AREA=0 

WRITE!  1,  1000  )ESC,AST 
FORMAT!  1X,2A1, 

' - -----JOYSTICK  CONTROL  MODES - " 

'  1.  MOVE  CURSOR  TO  SCREEN  CENTER'/ 

'  2.  MOVE  CURSOR  TO  PARTICLE  EDGE',' 

'  3.  AUTO  TRACK  AROUND  PARTICLE  EDGE'/ 

'  4.  DISPLAY  DIGITIZED  EDGE  OVER  PARTICLE'/ 

'  3+4.  RETURN  TO  OPERATING  SYSTEM'/ 

'  -X  POSN-Y  POSN— VIDEO-PTSIZE— -AREA— FERIft'/  ' 
GOTO  100 

WRITE!  l,1010)X0UT,Y0UT,VIDE0,PTSi:E,A!;EA,PtRI!l 
FORMAT! '  +  ',618) 

CALL  JSINP!JSK,JSX,JSY) 

IF!  JSK  '50,50,60 
IF!KJS)40.40.80 
IF!JSK-9)70,500,500 
KJS=JSK 

GOTO!  100, 200, 200, 300, 300, 300, 300, 400 ),  K  JS 
ZERO  CURSOR  OH  SCREEN 


X0UT=128 

Y0UT=128 

LOUT=130 

T0UT=0 

VINP=0 

AINP=0 

PINP=0 

CALL  TVSCN!XOUT,VINP) 

VIDEO=VINP 

Xl=128 

Yl=128 

6{I?5So 


MOVE  CURSOR  TO  PARTICLE  EDGE  WITH  JOYSTICK 


X1=X1+JSX/10 
Y1=Y1-JSY/10 
IF'Xl-254  )210.210. 30 
IF(Vl-254  '220,220,30 


1 


210 


220 

230 

240 


C 

C 

C 

300 

1300 

1310 

310 

1320 


320 

1330 


330 

1340 


340 


344 

346 


350 


360 


361 

362 

363 

364 


365 


370 

380 


1370 

1 


IF(Xl-2  ,'30.30,230 

IF<yi-2  >30,30, 240 

X0UT=X1 

Y0UT=Y1 

L0UT=Y0UT+2 

T0UT=0 

CALL  TVSCN(XOUT,VINP) 

VIDE0=VINP 

GOTO  30 

AUTO  TRACK  WAGE  EDGE 


MRITE(  1,1300) 

K'  ENTER  NUHBER  OF  SCAN  POINTS:  '  i 
,1310)NSCAN 
FORHATi 18 > 

IF(NSCAN-32)310,310,300 
WRITE(  1  1320) 

FORMAT! ENTER  SCAN  AMPLITUDE:  ') 

READ!  1, 1310  )SCAMF 
IF! SCAMP-16 >320, 320, 310 
WRITE!  1,1330) 

FORMAT! '  ENTER  SCAN  LEAD  START:  ' ) 

READ!  1, 1310  )SC  ART 

IF!  SCART-NSCAN/2  >330, 330, 320 

WRITE!  1,1340) 

FORMAT! ’  ENTER  UIDEQ  TOLERANCE:  ') 

READ!  1, 1310  )TOLER 

P0INT=NSCAN/4 

TW!0)=X0UT 

TH!0)=Y0UT 

TAREA=0 

TERIM=0 

JMP=0 

DA=6.28318/NSCAN 
DO  340  I=1,NSCAN 
DX!I)=SCAMP*COS!I»DA) 
DY!I)=SCAMP«SIN(I(!DA) 

INS=I+N5CAN 
DX!  INS)*[iX!I) 

DY!  INS  )=DY!  I ) 

CONTINUE 
DO  370  J=l,100 
START=POINT-SCART 
IF!  START  >344.346. 346 
START=NSCAN+P0INT-5CART 
DO  350  I=1.NSCAN 
IST=I*START 

XSCAN!l)sT«lJ-l)-DXlIST) 

YSCAN!I)=TH!J-1)+DY(IST) 

CONTINUE 


CALL  TVTRK! NSCAN, XSCAN!  1 ), YSCAN!  1 ) ) 

CALL  JSINP!  JSK,JSX,JSY) 
iF!JSK  >360,360,60 
TW!J)=XSCA/|!  POINT) 

TH!J>=YSCAN!  POINT) 

P0INT=P0INT+START 

IF! POINT-NSCAN >362,362, 361 

POINT=POINT-NSCAN 

IF!  3-2)370,365, 363 

IF!  IADS!  TW!  3  )-IW<  1 )  )tIA&S!  TH(  3  )-TH!  1 )  )-SCAMP+l  >364. 365, 365 
TW!3)=TW(1) 

TH!3)=TH!1) 


JMP'l 

DH=TH!3-1)-TH!3) 

DW=TW!3-1)-TW!3) 

AM=TW!3-l)+TW!3)-256 


TAREA=TAREA+AW»DH/2 
TERIM=TERIM+S!!RT!  Dll**2+(  FH»OH  mz  > 
IF!  3MP-1  >370,380,  330 
CONTINUE 
FTSIZE=3 

AREA=INT!FH*TAREA) 

PERIM*INT!  TERIM) 

WRITE! 1, 1370  )PTSIZE,PERIM. AREA 


FORMAT!'  NUMBER  PERIM  AREA  "318/ 
'  POINT  X  AXIS  Y  AXIS') 

DO  390  3M,PTSI2E 
WRITE!l.l380)3,TW!3),TH!J) 

FORMAT!  318) 


1380 


3?0 


CONTINUE 
GOTO  20 

C  DISPLAY  DIGITIZED  PERIMETER  OVER  IMAGE 

C 

400  T0UT=1 

DO  450  I=1,PTSIZE 

C  rnd=ran(i.o) 

C  J=INT(PTSIZE«RND) 

J=I 

XOUT=TU(  J) 

YOUT=TH(  J) 

LOUT=YOUT+2 

CALL  TVSCN(XOUT,VINP) 

CALL  JSINP(  JSK,JSX,JSY) 

IF(  JSK >450,450, 60 

450  CONTINUE 
GOTO  400 

C 

C  RETURN  TO  CALLING  PROGRAM 

C 

500  CONTINUE 

TL=PTSIZE 
TP=P£RIi1 
TA=AREA 
RETURN 
END 

A. 


TVTRK: 


SCOUT: 


STAT; 


NEXT: 


TOUT: 


TVTRK. Z80  SUBROUTINE 
REVISION  0.22,  4  AUGUST  19E2 

THIS  ROUTINE  MOVES  THE  CURSOR  TO  EACH  X  ANL  Y 
POSITION  AROUND  A  POINT  ON  THE  PARTICLE  EDGE 
UNTIL  A  CIRCLE  IS  COMPLETED  OR  THE  VIDEO  EXCEEDS 
THE  PREVIOUS  VIDEO  BY  A  CERTAIN  TOLERANCE. 

THE  POINT  AT  WHICH  THIS  OCCURS  IS  THEN  RETURNED. 

PARAMETER  ASSIGNMENT 

TVTRK(  NSCANCHLO, XSCANCDEO. YSCANCBC j ) 

YSCAN  ARRAY  HILL  START  AT  BC  AND  BE  RUT  IN  lY. 
XSCAN  ARRAY  HILL  START  AT  DE  AND  BE  RUT  IN  IX. 

NSCAN  ADDRESS  IS  IN  HL  REG.  VALUE  RUT  IN  E  REG. 
TOLER  IS  AT  HL+2,  VALUE  PUT  IN  C. 

EDGE  POINT  HILL  k  PUT  IN  HL+4  LOCATION  RGR  RETURN. 

SUBROUTINE  INITIALIZATION 


ENTRY 

TVTRK 

PUSH 

BC 

STORE  YSCAN  ADDRESS  IN  It 

POP 

lY 

PUSH 

DE 

STORE  XSCAN  ADDRESS  IN  IX 

POP 

IX 

LD 

E,(HL) 

HL 

STORE  NSCAN  VALUE  IN  £ 

INC 

INC 

HL 

LD 

C,(HL) 

STORE  VIDEO  TOLERANCE  In  C 

INC 

HL 

INC 

HL 

PUSH 

HL 

SAVE  POINT  ADDRESS  ON  STACK 

LD 

D,0 

INIT  POINT  IN  D 

LD 

6,0 

INIT  FOR  WHITE  TO  ElACK  EDGE 

OUTPUT  NEXT  SCAN  POINT 

INC 

D 

SETUP  FOR  NE<T  POINT 

LD 

A,  10 

LOAD  X  CONTROL 

OUT 

74H,A 

A, (IX) 

OUTPUT  X  CONTROL 

LD 

LOAD  X  VALUE 

our 

64H.A 

OUTRUr  X  VALUE 

LD 

A,  4 

LOAD  f  CONTROL 

OUT 

74H.A 

OUTPUT  I  CONTROL 

LD 

A,(IY) 

LOAD  t  value 

OUT 

64H,A 

OUTPUT  I  VALUE 

INC 

IX 

SETUP  POP  i)ExT  POSITION 

INC 

IX 

INC 

IT 

INC 

lY 

IN 

A,64H 

INPUT  STATUS 

BIT 

0,A 

CHECK  STATUS  BIT 

JP 

Z,STAT 

KEEP  CHECKING  UNTIL  STATUS  IS  GK 

IN 

A,74H 

INPUT  VIDEO 

CP 

B 

COMPARE  VIDEO  WITH  PREVIOUS  THRESHOLD 

JP 

C,POUT 

EDGE  IS  POUND  IP  VIDEO  IS  LESS 

CP 

C 

COMPARE  WITH  TOLER 

JP 

C,NEXT 

NEXT  POINT  IP  LESS 

SUB 

A,C 

SUB:TPACT  TOLER  FROM  Virso  IF  NOPE 

LD 

B,A 

AND  STORE  AS  NEXT  THPEShOLD 

DEC 

E 

DEC  NSCAN 

JP 

NZ, SCOUT 

REPEAT  IF  NOT  ZERO 

RETURN  EDGE  POINT 

LD 

A,D 

LOAD  EDGE  POINT  . 

POP 

HL 

LD 

m 

(HL).A 

AND  RETURN 

A. 


VIDEO  INTERFACE  3UPR0UTINE 
REVISION  0.21,  4  AUGUST  1932 

THIS  ROUTINE  MOVES  THE  CURSOR  TO  THE  CENTER  OF  THE 
TV  SCREEN.  PARTICLE  CENTER  OR  THE  PARTICLE  E^TRENES 
AND  SCANS  THE  PARTICLE  TO  FIND  THE  AREA  AND  FERINETEF. 

SUBROUTINE  PARAMETERS 

TVSCN(XINPCHL].V0UTCDE2) 

YINP  IS  AT  HL+2 
LINP  IS  AT  HL+4 
TINP  IS  AT  HL+6 
AOUT  IS  AT  DE+2 
POUT  IS  AT  DE+4 

DATA  BASE 


XINP: 

DB  1 

X  VALUE 

YINP: 

DB  1 

f  VALLE 

LINP: 

DB  1 

Y  LIMIT 

TINP: 

DB  1 

THRESHOLD 

VOUT: 

DB  1 

VIDEO 

AOUT: 

DB  1 

AREA 

POUT: 

DB  1 

PERIMETER 

ENTRY 

TV5CN 

TV3CN: 

LD 

A,(HL) 

BRING  IN  FIRST  PARAntTER 

ADD 

A.O 

SET  FLAGS 

JP 

Z.CNTR 

IF  ZERO  GO  CENTER  TV  CURE 

PUSH 

it 

SAVE  OUTPUT  ADDRESS 

LD 

A,(HL) 

STORE  X  VALUE 

LD 

(XINP),  A 

INC 

HL 

INC 

HL 

LD 

A,(HL) 

STORE  Y  VALUE 

LD 

(YINP),  A 

INC 

HL 

INC 

HL 

LD 

A,(HL) 

STORE  Y  LIMIT 

LD 

(LINP),  A 

INC 

HL 

INC 

HL 

LD 

A.  ( HL ) 

STORE  THRESHOLD 

LD 

(TINP),  A 

LD 

HL,CO 

RESET  HL 

LD 

DE.OO 

RESET  DE 

LD 

BC,00 

RESET  EC 

t 

OUTPUT  POSITION 

5iART: 

LD 

A,  10 

LOAD  X  CONTROL 

OUT 

74H,A 

OUTPUT  X  CONTROL 

LD 

A,  (XINP) 

LOAD  X  VALUE 
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ABSTRACT 


A  holographic  image  analysis  system  has  been  developed  to  measure  the 
position,  velocity,  size  and  shape  of  microscopic  particles  slowly 
settling  in  three-dimensional  space.  Images  of  particles  recorded 
sequentially  on  individual  holographic  frames  are  reconstructed  using  an 
in-line,  far  field  configuration.  Image  analysis  is  computer-controlled 
with  two  basic  functions.  First  is  the  precision  registration  and  XYZ 
positioning  of  the  holographic  frame  so  that  accurate  particle 
displacements  between  sequential  frames  can  be  measured.  These 
displacements  are  compared  to  the  elapsed  time  between  the  frames  to  get 
particle  settling  velocity.  Second,  it  scans  a  digitized  video  image  of 
the  reconstructed  holographic  image  to  measure  size,  shape,  and  area.  A 
scanning  algorithm  has  been  developed  to  determine  particle  size  (area, 
length,  width)  for  classification  by  settling  characteristics.  An  edge 
tracking  algorithm  has  also  been  developed  to  facilitate  particle 
identification  by  shape  from  frame  to  frame.  A  cataloging  system  was 


developed  to  provide  for  data  entry,  storage  and  retrieval. 
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I .  INTRODUCTION 

Measurement  of  size,  shape  and  density  are  critical  factors  in  the 
study  of  particle  dynamics  and  settling  characteristics  of  marine 
hydrosols.  Micrographic  holography  has  been  widely  used  in  aerosol 
studies  (1)  and  holographic  movie  cameras  have  been  developed  for  studying 
zooplankton  feeding  behavior  in  ocean  waters  (2).  The  size  and  settling 
behavior  of  aggregates  resulting  from  mixing  riverine  water  with  seawater 
has  been  holographically  determined  in  the  laboratory  (3)  and  effect  of 
particle  shape  on  the  settling  velocities  of  primary  hydrosols 
(unaggregated)  has  also  been  investigated  holographically  (4),  (5). 

Carder  e£  al^.  (6)  have  developed  an  ^  situ  holographic  mlcro- 
velocimeter  for  the  study  of  microscopic  particles  suspended  in  sea  water. 
This  device  records  the  image  of  each  particle  in  three-dimensional  space 
on  a  series  of  successive  holograms  with  respect  to  time.  Horizontal  and 
vertical  dimensions  and  cross  sectional  area  have  been  used  to  determine 
particle  density  using  Stokes'  Theorem  (3).  Edge  coordinates  are  useful 
for  particle  identification  (e.g.,  Zernike  moments)  and  analysis  of  shape 
on  particle  rotation  and  settling  speed  (4).  During  these  studies,  a 
great  number  of  samples  roust  be  analyzed  in  order  to  provide  sufficient 
statistics  for  determining  particle  population  characteristics  of  shape 
and  size  as  well  as  settling  velocity,  trajectory,  orientation  and 
oscillations.  Such  measurements  can  be  facilitated  by  use  of  a 


computer-controlled  image  analysis  system  to  Improve  the  accuracy  and 
reduce  the  number  of  man-hours  in  analyzing  holographic  particle  data. 


PROBLEM  DEFINITION 


II. 

A  measurement  technique  was  needed  to  determine  the  position,  size, 
shape,  orientation  and  velocity  of  microscopic  particles  moving  in  three- 
dimensional  space.  The  marine  particles  of  interest  in  the  study  (6) 
cited  above  typically  ranged  from  5  to  250  microns  in  diameter,  which 
necessitated  the  use  of  microscopic  techniques  for  any  detailed  size  or 
shape  analysis.  Since  oceanic  particle  densities  may  vary  from  about  1.03 
to  5.2  g/ml  and  typical  settling  velocities  may  range  from  less  than  .0001 
to  1.0  cm/s,  a  variety  of  frame  or  sample  periods  from  less  than  one 
second  to  as  much  as  one  minute  may  be  necessary  in  order  to  provide 
short-term  velocity  measurement  accuracy.  In  addition,  the  duration  or 
exposure  time  for  the  position  measurement  should  be  less  than  1/500  sec. 
in  order  to  maintain  sufficient  positional  accuracy  and  prevent  hologram 
"smearing."  Therefore,  a  relatively  high  speed,  microscopic  technique  was 
developed  (6)  to  record  the  images  of  multiple  particles  as  they  settled 
in  three-dimensional  space.  In  order  to  deal  with  the  enormous  data 
volume  generated  in  applying  these  techniques  (more  than  50  particles  per 
hologram),  an  automated  method  for  reconstructing  and  analyzing  the 
particle  data  has  been  developed.  Edge  coordinates  are  recorded  when 
needed  for  particle  identification  and  for  analysis  of  the  particle  shape 
effects  on  particle  rotation. 

III.  SYSTEM  DESCRIPTION 

A  functional  block  diagram  of  the  holographic  image  analysis  system 
is  showm  in  Figure  1.  The  major  components  and  functions  are:  1)  a  laser 
to  reconstruct  the  holographic  image,  2)  a  mechanical  stage  to  position 
the  hologram  in  three-dimensional  space,  3)  a  magnification  lens  and  rear 
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i 


2)  XYZ 


FIGURE  1  HOLOGRAPHIC  IMAGE  ANALYSIS  SYSTEM 
BLOCK  DIAGRAM 


projection  screen  on  which  to  reconstruct  the  particle  image,  4)  a 
television  camera  to  view  the  images,  5)  a  video  digitizer  to  convert  the 
image  into  digital  data,  6)  a  television  monitor  to  display  the  particle 
image,  7)  a  microcomputer  with  disk  storage  for  processing  and  storing 
input  data  from  the  video  digitizer  and  controlling  the  position  of  the 
mechanical  stage,  8)  an  interactive  terminal  to  enable  operator  control  of 
the  system,  9)  a  joystick  for  manual  control  of  the  frame  position  and  10) 
the  XYZ  drive  controller  to  move  the  holographic  frame  in  the  laser  beam. 
Table  1  presents  a  list  of  the  equipment  which  comprises  the  system. 


Table 

1.  HOLOGRAPHIC  IMAGE  ANALYSIS  SYSTEM  EQUIPMENT 

I  tern 

Description 

Manufacturer 

Model 

I 

Laser,  15  mK,  He-Ne 

Spectraphysics 

124B 

2 

XYZ  axis  table 

Aerotech,  XY : 

ATS-303 

Z: 

ATS-416 

3 

28rim  or  50mm  lens 

Nikon 

4 

Television  camera,  high  resolution 

Dage 

650 

5 

Video  digitizer 

Colorado  Video 

270A 

6 

Television  monitor 

Panasonic 

W  5300 

7 

Micro-computer,  Z-80,  dual  disk 

Cromemco 

Z-2D 

8 

CRT  terminal 

Soroc 

IQ  120 

9 

Joystick 

Cromemco 

10 

XYZ  drive  controller 

Aerotech 

EC-2 

The  system  performs  the  following  functions:  1)  determines  the 
position  of  a  particle  with  respect  to  a  three-dimensional  frame 
reference,  2)  determines  the  width,  height,  cross  sectional  area  and 
perimeter  of  a  particle  using  a  scanning  algorithm,  3)  determines  the 
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cross-sectional  area,  perimeter  and  edge  coordinates  of  the  particle  using 
a  tracking  algorithm,  4)  provides  for  particle  identification  from  frame- 
to-frame  by  displaying  the  particle  position  and  edge  points  (shape)  from 
the  previous  holographic  frame  on  the  current  frame  and  5)  computes  the 
velocity,  average  diameter  and  particle  density  from  the  recorded  data  for 
each  particle. 

IV.  TECHNICAL  DISCUSSION 

Successive  holograms  of  settling  particles  were  recorded  iai  situ  on 
negative  holographic  transparency  film  using  the  particle  velocimeter 
described  in  reference  (b).  The  film  was  processed  using  standard 
photographic  techniques  and  mounted  on  individual  frames.  Each  holo¬ 
graphic  frame  was  positioned  In  a  laser  beam  to  reconstruct  a  two- 
dimensional  profile  of  each  particle  in  the  direct  transmission  mode.  The 
image  of  each  particle  was  successively  focused  on  the  rear  projection 
screen  by  moving  the  holographic  frame  in  the  Z  axis.  The  projected  scene 
was  viewed  with  a  television  camera  and  converted  into  a  digitized  video 
signal.  The  holographic  frame  was  then  moved  in  the  X  and  Y  axes  to  bring 
the  focused  image  of  each  particle  into  the  center  of  a  TV  monitor.  The 
position  of  the  frame  and  the  size,  shape  and  orientation  of  each  particle 
was  analyzed  and  recorded  by  the  computer.  As  images  of  the  sarnie  parti¬ 
cles  in  subsequent  frames  were  catalogued,  a  matrix  of  settling  velocity 
and  lateral  motion  characteristics  was  generated  for  each  particle.  The 
statistical  distribution  of  settling  velocity  with  respect  to  particle 
size,  shape  ,'”id  orientation  provided  an  accurate  determination  of  particle 
density  and  ultimately  permitted  a  general  particle  classification  (e.g., 
organic,  mineral,  heavy  mineral),  based  largely  upon  density. 
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A. 


PARTICLE  POSITION  AND  VELOCITY 


Particle  settling  velocity  using  the  system  described  in  Figure  1  was 
determined  by  measuring  the  position  of  each  particle  in  three-dimensional 
space  at  two  or  more  different  sample  times.  A  typical  particle 
configuration  is  shown  in  Figure  2.  Prior  to  cataloging  any  particle 
images  from  a  frame,  an  accurate  and  repeatable  reference  is  established 
for  a  point  common  to  all  frames.  The  reference  point  in  a  frame  is 
centered  on  the  television  monitor  by  moving  the  frame  in  the  X,  Y  and  Z 
axes  under  computer  control  with  a  joystick.  The  image  of  each  particle 
in  the  frame  is  then  subsequently  moved  to  the  center  of  the  monitor  with 
the  Joystick  and  its  position  recorded.  Size  and  orientation  parameters 
are  determined  using  the  particle  scanning  and  edge  tracking  algorithms 
described  below.  Once  recorded,  the  particle  data  from  a  previous  frame 
can  be  retrieved  and  the  XYZ  stage  driven  to  that  position  under  computer 
control  to  enhance  location  and  identification  of  that  particle  or.  a 
subsequent  frame.  When  the  entry  of  data  for  all  frames  has  been 
completed,  the  distance  each  particle  has  travelled  between  successive 
frames  is  determined  and  the  individual  velocities  calculated  as  shown: 

Vt  =  (Vx^  +  Vy^  +  (1) 


where 

Vx  =  dx/dt 
Vy  =  dy/dt 
Vz  =  dz/dt 
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B.  PARTICLE  AREA  AND  PERIMETER 


DeCermination  of  the  projected  cross-sectional  area  of  a  particle  can 
be  derived  from  the  reconstructed  image  by  counting  the  number  of  pixels 
inside  the  particle  boundary.  In  a  similar  manner,  the  perimeter  of  the 
particle  can  be  derived  from  the  number  of  pixels  along  the  edge.  One  of 
the  major  problems  with  a  "noisy"  image  is  detection  of  the  pixels  which 
represent  the  edge.  Yakimovsky  discussed  an  edge  detection  technique  (7) 
based  on  maximizing  the  likelihood  ratio  with  a  simple  single-pass, 
region-growing  algorithm.  This  ratio  is  a  comparison  of  the  Intensities 
of  neighboring  pixels  within  the  same  object  to  neighboring  pixels  from 
two  different  objects.  The  boundary  derision  is  based  on  comparing  the 
ratio  with  a  pre-defined  threshold  and  accumulating  those  neighborhoods 
which  are  considered  to  be  within  a  single  object  into  one  contiguous 
region.  McKee  and  Aggarwal  have  developed  a  multi-pass  technique  (8) 
which  processes  a  full  video  frame  of  data  to  obtain  edge  coordinates  of 
very  complex  shapes.  However,  this  method  would  be  far  more  complex  and 
time-consviming  than  necessary  for  determination  of  area  and  perimeter  in 
most  of  our  applications.  Also,  it  is  rare  for  more  than  one  particle 
image  to  be  found  in  a  given  focal  plane. 

Therefore,  we  have  chosen  to  use  a  simplified  edge  detection 
technique  based  upon  a  comparison  of  the  pixel  intensity  with  a  threshold 
adjusted  for  the  average  intensity  between  the  particle  and  background 
(SCAN  technique)  and  a  simplified  edge-tracking  method  (TRACK  technique)  . 
Pixels  at  each  transition  across  the  threshold  are  accumulated  as  the  edge 
of  the  particle,  while  pixels  above  the  threshold  are  accumulated  as  its 
area . 
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C.  PARTICLE  AREA  SCAN  TECHNIQUE 

The  particle  scanning  algorithm  developed  for  this  system  is  based  on 
a  simplified  single-pass  edge  detection  process.  The  surface  area  and 
edge  perimeter  of  each  particle  are  determined  by  comparing  the  video 
intensity  of  each  pixel  with  a  selected  threshold  value.  The  technique 
consists  of  processing  a  series  of  vertical  scans  which  horizontally  cover 
the  image  as  depicted  in  Figure  3.  The  horizontal  (equatorial)  and 
vertical  (polar)  axes  of  the  particle  are  determined  by  manually  adjusting 
a  vertically  scanning  cursor  to  the  edges  of  the  particle  when  the 

algorithm  is  first  executed  for  each  measurement.  This  limits  the  region 
under  investigation  to  the  immediate  environment  of  the  particle  and 
eliminates  the  generally  out-of-focus  neighboring  particles  from 
consideration . 

Beginning  on  the  loft  side,  the  particle  is  scanned  from  top  to 

bottom  in  synchronism  with  the  television  raster.  The  video  intensity 

(0-255)  of  each  pixel  is  compared  with  a  threshold  value  corresponding  to 

the  intensity  at  the  particle  edge.  If  the  pixel  intensity  is  greater 

than  the  threshold,  the  pixel  area  is  assigned  a  unit  value  and  the  column 

area  is  incremented.  Otherwise,  its  value  is  zero.  If  the  pixel  area 

value  is  different  from  that  of  the  last  pixel,  the  horizontal  perimeter 

is  incremented.  When  the  vertical  scan  is  complete,  the  difference 

between  areas  of  the  ctirrent  and  previous  columns  is  added  to  the  vertical 

perimeter.  The  scan  is  then  incremented  to  the  right  and  the  process 

repeated  until  the  entire  particle  has  been  scanned.  The  selected  scan 

consists  of  N  rows  vertically  and  N  columns  horizontally.  The  total 
r  c 

area  of  the  particle  is  the  sum  of  all  column  areas  scaled  as  shown  in  Eq . 
(2).  Since  the  shape  of  each  pixel  is  not  square,  the  vertical  component 
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must  be  scaled  to  provide  an  accurate  measurement  of  the  area  and 
perimeter.  Or  simply. 

If 


If 


V  / .  .  ^  then  A  , .  . .  =  I 

p(i.j)  t  p(i,j) 


V  V  then  A  ..  ..  =  0 

P(i.j)  t  P(i.j) 


N 


N 


(2) 


where 


Area  =  F  *  F  *  Sum  (  Sum  A  (i,j)  ) 

X  V  .  .  1  P 

1=1  J=1 

i  and  j  are  the  pixel  column  and  row, 

V  is  Che  video  amplitude  of  pixel  (i,j), 

F 

is  Che  preselected  video  threshold, 

N  is  the  number  of  seamed  coliim.ns, 
c 

is  the  number  of  scanned  rows, 

F  is  the  horizontal  scale  factor  and 

X 

F  is  the  vertical  scale  factor. 


The  error  introduced  by  the  Quantization  of  area  should  converge  toward 
zero  and  become  negligible  as  the  size  of  the  particle  increases. 

The  vertical  component  of  the  perimeter  »  is  the  difference 

between  Che  number  of  pixels  above  the  threshold  in  one  column  compared  to 
the  number  in  the  previous  column.  The  horizontal  component  (Pj,) >  i® 
Che  sum  of  all  transitions  across  the  threshold  occurring  in  each  column. 
The  horizontal  and  vertical  components  are  scaled  by  F^  and  F^, 
respectively,  and  combined  to  form  the  particle  perimeter  as  shown  in 
Eq.  (3).  If  Vp(i,j)  is  not  equal  to  V^(i,j  +  1)  then  the  Incremental 
horizontal  component  (il  (i,i)  =  I,  otherwise  H  (i.i)  =  0. 

p  -  p  ” 


Perimeter  =(F  *P.  )+(F  *P  ) 

X  h  V  V 


(3) 
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whore 


and 


Ph 


K 

-  Sum  ( 
i=*l 


t; 

Sum  H  (i,i) 

j  =  l  P  ' 


N  N  N 

P  *  Sum  (  Sum  A  (i,i)  -  Sum  A  ). 

i=l  j  =  l  P  j  =  l  P 
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The  major  error  in  this  method  is  highly  dependent  on  the  orientation  of 
any  edge  segment  with  respect  to  the  scanning  axis.  Although  straight 
edges  orthogonal  to  the  scan  will  yield  good  measurements  of  perimeter, 
those  sane  straight  edges  oriented  diagonally  to  the  scan  axis  can  degrade 
the  perimeter  measurement  by  as  much  as  30%  [1-SQRT  (1/2)].  Therefore,  a 
different  method  of  determining  perimeter  which  will  be  Insensitive  to 
particle  orientation  was  required.  Also,  variations  in  the  pixel 
intensity  around  the  particle  perimeter  with  respect  to  the  threshold 
level  increase  measurement  uncertainty.  A  dynamic  threshold  which  can 
adapt  to  varying  particle  and  background  video  intensities  will  improve 
this  condition.  Such  a  technique  is  implemented  in  an  edge-tracking 
routine  described  below  that  provides  data  for  determination  of  particle 
shape  and  orientation. 


D.  PARTICLE  SHAPE  AND  ORIENTATION 

Boundary  definition  using  direction  and  curvature  chains  have  been 
described  by  Eccles  e^  aj^.  (9).  These  techniques  generally  describe  an 
image  boundary  in  terms  of  a  one-dimensional  list  of  angles  starting  from 
an  origin  on  the  Image  edge  and  incrementing  uniformly  around  the 
boundary.  Freeman  has  described  another  method  for  determining  shape  by 
defining  critical  points  around  the  boundary  (10).  These  points  Include 
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discontinuities  in  curvature,  points  of  inflection,  curvature  maxima, 
intersections  and  points  of  tangency.  Fourier  transforms  have  been  used 
as  digital  filters  to  smooth  the  digitized  boundaries  of  planar  objects. 
Pavlidis  discusses  several  algorithms  using  Fourier  transform  coefficients 
to  define  detailed  shape  characteristics  (11).  Teague  has  developed 
similar  shape  analysis  techniques  using  Zernike  moments  (12). 

We  have  chosen  a  boundary  chain  technique  to  determine  particle  shape 
and  orientation  which  is  similar  to  the  method  described  by  Eccles  (9). 
The  edge-tracking  algorithm,  described  below,  determines  the  boundary 
between  the  particle  and  background  by  examining  a  small  region  of  pixels 
along  the  particle  edge.  Edge  detection  is  based  on  the  intensity 
gradient  between  neighboring  pixels  rather  than  the  simple  threshold 
technique  used  in  the  particle  scanning  algorithm  described  above.  As 
each  new  point  along  the  the  particle  edge  is  determined,  the  tracking 


pattern  ii 

ri  then 

centered  on  that 

point 

and 

the  new 

neighborhood  is 

explored . 

This 

process  continues 

until 

the 

particle 

has  been  fully 

circumnavigated.  A  more  detailed  discussion  is  given  below. 

E.  PARTICLE  EDGE-TR.\CKING  TECHNIQUE 

The  particle  edge-tracking  algorithm  provides  the  capability  to 
classify  particle  size,  shape  and  orientation  for  applications  dealing 
with  particle  settling  dynamics.  The  edge  tracking  algorithm  scans  across 
the  particle  edge  in  a  series  of  small  circular  patterns  Xv’here  i  is  one  of 
the  Np  scans  around  the  particle  edge  and  j  Is  one  of  the  N^  points  in 
a  scan  as  shown  in  Figure  4.  The  point  j  in  the  scan  for  which  the  pixel 
Intensity  V(i,j)  exceeds  the  previous  pixel  intensity  V(l,j-1)  by  a 
threshold  value  V^,  is  defined  as  the  edge  for  the  current  scan  pattern 
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i  and  is  used  as  the  center  of  the  next  scan  pattern  (i+1).  This  process 
continues  until  the  scan  returns  to  the  neighborhood  of  the  starting 


point.  The  scan  radius  R^,  number  of  points  per  scan  and  threshold 
value  V^,  can  be  manually  selected,  depending  upon  the  particle  size  and 
shape. 

The  perineter  of  each  particle  is  measured  by  accumulating  the  radial 
distance  of  each  threshold  transition  around  the  particle  edge  until  the 
circumnavigation  of  the  particle  is  complete.  Since  the  shape  of  each 
pixel  is  not  square,  both  vertical  and  horizontal  components  must  be 
scaled  accordingly,  by  and  F^,  respectively,  to  provide  an  accurate 
measurement  of  the  distance  around  the  particle.  The  algorithm  for  this 
function  is: 


N 

P  I 

Perimeter  =  Sum  (  (  F  *  T  )==  +  (F  *  T  ]  '"  (k) 

.  ,  y  y  X  X  ^  ' 

1=1 

where  T  =  R  sin 
y  s 

T  =  R  cos 
X  s 

R  =  scan  radius 
s 

Np  =  number  of  threshold  transition  points  along  edge. 

The  area  of  each  particle  is  determined  by  accumulating  the 
trapezoidal  area  defined  by  eacli  new  transition  as  shown  in  Figure  A.  The 
distance  between  the  particle  edge  and  a  vertical  reference  line  at 
the  center  of  the  screen  is  determined  for  each  transition  point  and 
multiplied  by  the  vertical  distance  dy^  between  two  successive  edges. 
The  entire  area  of  the  particle  is  determined  by  successively  adding  or 


subtracting  each  incremental  area  depending  on  the  vertical  direction 
until  the  particle  has  been  completely  encircled  as  described  by: 


N 

P 

Area  =  (  1/2  Sum  (x+x  )*dy*F/F)  (5) 

i=l  i  i+1  i  ^  ^ 

Particle  identification  from  frame-to-f rame  is  achieved  by  displaying  the 
edge  pixels  of  a  selected  particle  on  the  next  frame  and  comparing  it  to 
the  shapes  of  the  new  particles  for  correlation. 

E.  ALGORITHM  ERROR  ANALYSIS 

The  inherent  error  in  the  tracking  process  consists  of  radial  and 

angular  coinponents.  as  shown  in  Figures  5  and  6,  respectively.  These 

errors  are  functions  of  the  number  of  scan  points  »  '"’tid  the  scan 

radius  (R  ),  as  shown  in: 
s 

Radial  error  =  R  *  (  1  -  cos  (2  -r  /  N  ))  /  2 
s  s 

Angular  error  ==  +/-  2  '  /  Ng 

These  errors  on  average  will  produce  a  small  positive  error  (e.g., 
2.5”<  for  N  =  20)  in  determining  the  perimeter  of  a  smooth  boundary  as 
shown  in  Eq.  (7).  However,  the  algorithm  may  also  exhibit  a  negative 
"corner  truncation"  error  for  figures  with  sharp  corners  because  of 
truncation  produced  by  the  discrete  scan  process.  The  difference  between 
the  actual  and  measured  area  Cor  large,  smooth  particles  should  be  small 
but  may  become  negative  for  particles  with  sharp  corners,  Eq.  (8). 
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FIGURE  5  RADIAL  ERROR  DUE  TO  AS 
AFFECTING  PERIMETER  ACCURACY 


FIGURE  6  ANGULAR  ERROR  AFFECTING 
AREA  ACCURACY 


(7) 


Peritn  error  =  +  N  *R  (1-  cos  2  tt  /  N  ) 
p  s  s 

Area  error  =  +/-  (  N  *  R  *  sin  2  /  N  )  /  2  (8) 

p  s  s 

Area  and  perimeter  measurements  will  also  vary  as  a  function  of  the 
video  intensity  threshold  value  for  the  scan  mode  and  as  a  function  of  the 
video  intensity  gradient  threshold  value  for  the  tracking  mode.  The  video 
amplitude  across  the  edge  of  an  image  with  varying  contrasts  is  shown  in 
Figure  7.  The  effects  due  to  variation  of  the  threshold  values  on 
measurements  of  a  high-contrast,  white-on-black  square  (2  cm/side)  are 
summarized  in  Table  2.  In  generating  these  data  the  video  threshold 
intensity  for  the  scar,  mode  was  set  at  values  from  237,  to  IWA  of  the 
intensity  difference  between  the  square  and  the  background,  whereas  the 
intensity  gradient  chre.shold  (intensity  difference/pixel)  for  the  tracking 
algorithm  was  varied  from  87;  to  327(  of  the  total  Intensity  difference 
between  Che  square  and  the  background.  This  particular  application  used  a 
negative  threshold  gradient  (i.e.,  a  bright-to-dark  transition). 

Tabic  2.  Effects  of  variation  of  video  threshold  intensity  (scan)  and 

intensity  gradient  (track).  Variation  is  represented  by  ±  one 
standard  deviation  as  well  as  the  total  range  of  deviation. 

AREA  (PIXELS)  PERIMETER  (PIXELS) 

Scan  Track  Scan  Track 

2165.6  ±  4.0%  2058.1  ±  1.3%  185.6  +  4.7%  180,4  i  1.77; 

llZ  Range  4%  Range  157o  Range  5.57(  Range 

The  range  of  threshold  effects  is  smaller  for  the  traclclng  mode  than 
the  scanning  mode  because  the  maximum  video  intensity  gradient,  generally 
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VIDEO  AMPLITUDE  VIDEO  AMPLITUDE 


SELECTED  PIXEL  =  178  OR  204 


HIGH 

CONTRAST 


BACKGROUND  PARTICLE  PARTICLE 
EDGE  INTERIOR 


SELECTED  PIXEL  =  118  OR  82 


FIGURE  7  VIDEO  AMPLITUDE  ACROSS  THE  EDGE  OF  A  HIGH 
CONTRAST  (TOP)  AND  LOW  CONTRAST  (BOTTOM) 
PARTICLE 


found  about  halfway  between  the  intensities  of  the  image  and  the  back¬ 
ground  (see  Fig.  7),  cannot  be  exceeded  without  the  algorithm  "losing 
track"  of  the  particle  edge.  The  threshold  intensity  gradient  value  for 
the  tracking  program  is  generally  selected  to  be  smaller  than  the  smallest 
normal  edge  gradient  (bottom  curve  in  Fig.  7)  found  along  the  particle 
edge.  This  insures  that  oblique  approaches  to  the  edge  do  not  reduce  the 
measured  intensity  gradient  below  that  of  the  threshold.  This  reduced 
threshold  value  causes  the  tracking  program  for  the  high-contrast  part  of 
the  particle  (upper  curve.  Figure  7)  to  in  effect  underestimate  the  size 
of  the  particle  for  a  white-on-black  image  and  overestimate  the  size  of  a 
particle  for  a  black-on-white  image,  especially  for  particles  of  varying 
contrast  along  the  edges.  For  the  white  square  on  black  background  image 
addressed  in  Table  2,  this  resulted  in  particle  size  estimates  that  were 
h  to  Ih  pixels  lower  per  edge  transition  (see  upper  curve  in  Figure  7)  , 
depending  upon  the  video  threshold  intensity  gradient  selected,  or  a  1.2 
to  3.6%  underestimate  of  the  length  of  each  side.  This  underestimate 
would  be  offset  to  some  extent  by  the  radial  error  effect  on  the  perimeter 
as  described  by  Eq.  (7). 

IV.  PERFORMANCE  EVALUATION 

The  performance  of  the  scan  and  track  algorithms  was  evaluated  using 
white-on-black  geometrical  figures  (square,  equilateral  triangle,  and 
circle)  of  known  size.  The  standards  for  comparison  were  accurate  to 
within  h  to  1  pixel  per  side  or  to  about  ±2.5%  for  the  small  figures  and 
±1.25%  for  the  large.  The  results  of  the  comparisons  are  shown  in  Table 
3.  Digitization  error  appears  to  have  produced  a  greater  effect  on 
accuracy  of  the  small  figures  with  complex  shapes  compared  to  the  larger 
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ones.  A  digitization  error  of  one  pixel  per  edge  can  result  in  errors  as 
large  as  ±3.6%  for  the  small  and  ±2.0%  for  the  large  triangle 
measurements.  However,  this  effect  should  approach  zero  as  the  shapes 
become  larger  and  smoother. 


Table  3.  Perimeter  and  area  measurement  accuracies  of  the  scanning  and 
tracking  algorithms  for  some  standard  geometrical  shapes 

Track  Scan 


Shape 

Size 

(cm) 

Area 

%  Error 

Perimeter 

%  Error  1 

Area 

:  Error 

Perimeter 

%  Error 

Square 

3.0 

0.38 

-0.47 

+0.50 

+  1.4 

Square 

2.0 

0.98 

0.17 

+  2.89 

+3.2 

Triangle 

4.0 

-3.00 

-2.90 

-0.56 

6.89 

Triangle 

2.2 

+3.98 

+  1.83 

-0.80 

+4.23 

Circle 

4.0 

-0.90 

0.62 

+  3.14 

+14.45 

Circle 

2.2 

-0.88 

2.97 

+0.66 

+9.36 

% 

Difference 

For  the  tracking  program 

,  most  errors  did  not  exceed 

the  accuracy 

limits  of 

the  standard  measurement.  Digitization 

errors, 

together  with 

the  uncertainties 

in  the  standards,  can  account 

for  the  differences 

between  the  measured  and  "standard"  sizes;  inclusion  of  the  potential  of 
radial,  corner-truncation,  and  threshold  intensity/threshold  intensity 
gradient  errors  insures  that  the  maximum  tracking  error  measured  (3.98%) 
falls  within  potential  error  limits. 

For  the  scan  program,  all  of  the  areas  measured  were  within  expected 
error  limits,  although  the  triangle  and  circle  perimeter  errors  were 
larger.  The  perimeter  errors  for  the  squares  were  quite  low  as  the  edges 


15 


of  each  square  were  aligned  with  the  vertical  and  horizontal  axes  of  the 
TV  camera.  Perimeters  in  the  scan  mode  can  be  exaggerated  by  as  much  as 
1  -  ''2/2  for  lines  sloping  45“  to  the  vertical.  The  largest  manifestation 
of  this  "serrated  edge"  effect  was  found  for  the  circle,  with  more  than  a 
+14%  error. 


V.  APPLICATION  TO  HOLOGRAMS 

To  demonstrate  the  application  of  the  system  and  techniques  to  the 
analysis  of  situ  transmission  holograms,  they  were  utilized  to 
re-examine  holographic  images  of  mica  flakes  collected  during  a  settling 
experiment  reported  in  (5) .  Because  the  mica  flakes  are  very  thin  and 
tend  to  settle  at  some  angle  to  the  focal  plane,  only  one  side  will  be  in 
sharp  focus  at  a  time.  Therefore,  a  slightly  defocused  (lower  contrast) 
image  was  analyzed  (Figure  8).  This  lower  contrast  image  is  typical  of 
the  type  collected  from  irregularly-shaped  particles  commonly  found  in 
seawater. 

A  comparison  between  the  tracking  and  scanning  algorithm  data  (Table 
4)  shows  that  the  areas  generated  are  within  about  6%  of  each  other  for 
these  mica  holograms.  The  perimeter  scanning  algorithm  serves  best  as  an 
upward  limit  or  check  on  the  perimeter  tracking  algorithm.  For  example, 
the  circular  image  of  the  third  holographically  recorded  particle  resulted 
in  a  perimeter  that  was  11%  higher  for  the  scanning  than  for  the  tracking 
algorithm.  Since  this  was  a  circular  particle,  the  scanning  algorithm 
overestimated  the  perimeter  by  serrating  the  edge. 
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Figure  8. 


Photograph  showint;  video  displav  of 
holoi;ram.  Two  ri.;ht-mosL  bars  are 
scale  for  all  pixels  alon-  central 


reconstructed  mica 
0-256  video  amplitiui 
bar . 


f  lake 


Table  4.  Holographic  inage  analysis  measurements  of  falling  mica  flakes 
using  the  scanning  and  tracking  algorithms  (magnification  of 
150X) 

- Scanning -  - Tracking - 

Threshold  General 


Sample 

Area 

Perlm 

Thresh 

Area 

Perim 

Gradient 

Shape 

1 

498 

127 

134 

500 

135 

26 

Rectangular 

2 

107 

40 

113 

101 

42 

30 

Rectangular 

3 

485 

118 

100 

501 

106 

30 

Circular 

VI .  SUMMARY 

The  holographic  microvelocimeter  has  proven  to  be  an  invaluable  tool 
in  studying  the  settling  dynamics  of  particles  both  in  the  marine 
environment  and  in  the  laboratory.  The  data  reduction  of  this  tool  has 
been  semi-automated  using  a  microcomputer  to  provide  data  bs.se  management 
and  to  control  two  phases  of  the  analysis.  In  the  first,  a  set  of 
precision  translation  stages  register  different  holographic  frames  and 
measure  particle  displacement  between  frames.  The  accuracies  of  this 
stage  are  on  the  order  of  ±0.0002  cm.  Even  greater  accuracies  in  particle 
settling  velocities  can  be  achieved  by  adjusting  the  exposure  interval. 

The  second  microcomputer-controlled  phase  is  the  measurement  of  the 
size,  shape  and  perimeter  of  the  digiti.zed  video  image  of  the  hologram. 
Two  techniques  have  been  developed  to  measure  the  area  and  perimeter  of  an 
image.  The  first  is  a  scanning  algorithm  which  utilizes  an  absolute 
threshold  of  the  video  intensity  of  each  pixel  inside  the  video  image. 
This  is  an  extremely  fast  algorithm  which  is  able  to  scan  the  area  and 
perimeter  in  a  single  pass.  It  suffers  in  accuracy  when  dealing  with 
lower  contrast  particles  having  a  variable  Intensity  in  different 
sections,  and  witli  rounded  particles  since  it  tends  to  count  a  serrated 
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edge.  A  second  technique  has  been  developed  to  track  the  edge  of  an  image 
and  increment  the  area  and  perimeter  utilising  an  edge  gradient  between 
the  background  and  image.  Both  algorithms  have  been  tested  with  standard 
geometric  images  and  with  holographic  images  and  are  accurate  to  within 
4-6%. 

While  this  system  has  been  developed  to  study  particle  settling 

dynamics,  the  techniques  of  image  analysis  could  be  applied  to  many 
pattern  recognition  and  any  digital  image  analysis  problems.  The  edge 
coordinates  from  the  cracking  algorithm  can  be  saved  and  further  analysis 
of  shape  done  using  Fourier  and  Zernicke  moments. 

This  effort  was  funded  under  Office  of  Naval  Research  contract 

NOOO 14-75-C-0539  to  the  Marine  Science  Department,  University  of  South 
Florida. 
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ABS'l'RACT 


SeCLlinc;  experiments  performed  on  silt  to  fine  sand  sized  mica  flases 
with  a  ho  lop.rapliic  micro-velocimeter  revealed  that  mica  is  the  hydraulic 
equi\’alent  of  quartz  spheres  having  diameters  a  factor  of  4  to  12  times 
siitalier.  Mica  in  the  very  fine  to  fine  sand  sizes  has  been  traditionally 
used  by  sedimentoiogists  to  delineate  areas  of  deposition  or  non-deposi- 
tion  and  potential  winnowing  of  fines,  and  is  here  found  to  be  the  hydraulic 
equivalent  of  silt  sized  particles  but  not  of  clay. 

Experiments  also  showed  that  mica  flakes  tend  to  settle  at  orientations 
which  are  neither  perpendicular  nor  parallel  to  the  gravitational  vector 
and  to  generally  maintain  their  orientation  throughout.  Equations  for  the 
settling  of  a  disc  in  Lerman  and  others  (1974)  and  that  developed  by  Komar 
and  Reimers  (1973)  are  shown  to  be  mathematically  similar  for  the  coarse 
silt  to  fine  sand  ranges  of  discs  and  are  adequate  predictors  of  settling 
rates  of  mica  flakes.  A  comparison  of  the  hydraulic  equivalency  of  quartz 
solieres  to  coarse  silt  througli  fine  sand  sized  mica  flakes  is  nresenfed. 
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INTUODUC!  ION 


liocausc  it  cleaves  into  flaky  partieJes,  sand  sized  mica  is  a  mineral 
i;roup  in  vhicii  shape  should  obviously  affect  settling  chnr.ic  t  e  r  is  t  ic  s  and, 
in  tact,  has  ol ten  been  considered  to  be  the  hydraulic  equivalent  of  silt 
ami  clay.  N'eiliesel  (1965)  noted  die  close  association  between  sand  sized 
mica  and  tiie  clay  fraction  of  Georgia  estuaries,  an  association  also  recog¬ 
nized  by  Panerancblum  (1966)  off  Israel.  Uoyle  and  others  (1968)  used 
the  abundance  of  mica  in  the  125-250  pm  size  fraction  to  delineate  areas 
of  the  southeastern  United  States  continental  margin  which  might  be  under¬ 
going  winnowing  or  deposition  of  fines,  a  process  which  otherwise  would  be 
masked  by  the  dominance  of  reworked  Pleistocene  sands.  A  similar  approach 
was  used  by  Adegobe  and  Stanley  (1972)  on  the  Niger  Shelf.  Doyle  and  othe 
(1979)  and  Park  and  Pilkey  (1981)  discuss  the  significance  of  mica  content 
to  the  deposit ional  and  erosional  systems  of  the  whole  continental  margin 
of  the  Eastern  United  States. 

Despite  its  intuitive  widespread  use  as  a  hydraulic  analog  of  finer 
sized  sediments,  no  quantitative  evaluation  of  the  sedimentologic  charac¬ 
teristics  of  mica  has  yet  been  undiirtakon.  The  purpose  of  this  paper  is 
to  determine  as  far  as  possible  the  liydraulic  characteristics  of  mica 
flakes  and  how  tliese  characteristics  compare  with  those  predicted  from 
settling  equations  in  the  literature. 

APPKOACIl 


Our  approacli  is  to  utilize  a  modifiei)  holograohic  micro-velociineter 
developed  and  described  hy  C.irder  (1978)  and  Carder  and  Meyers  (1979). 
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I'i.-.uio  !  siiuu’s  Lliu  sy.stoni  used.  Transmission  holograms  are  simultan- 
eoiislv  i-ollected  along  the  vertical  and  horizontal  axes  of  a  settling 
cuvi-itc.  A  reference  point  in  the  cuvette  is  also  recorded  on  each 
vortical  and  horizontal  hologram  to  allow  translation  between  the  two 
axes.  Seiiuential  frames  at  accurately  timed  intervals  (from  0.5  to 
J.i)  secoiuls  depending  on  the  particle  size)  record  the  settling  velocity 
■itui  orieni.iLion  as  well  as  particle  size  and  shape.  The  images  are 
reconstructed  by  placing  the  hologram  back  into  the  laser  path  and  re¬ 
focusing  on  particles  in  any  of  the  infinite  number  of  focal  nlanes  in  the 
settling  cuvette.  Tliis  system  offers  the  advantage  of  using  actual 
sedimentary  particles  in  the  sand  to  clay  size  ranges,  thus  obviating 
any  errors  which  may  be  inherent  in  scale  modeling  systems. 

.'lica  samples  were  chosen  from  sediment  cores  from  the  Eastern  United 
St.Ues  continental  slope  (Doyle  and  others,  1979).  Silt  to  fine  sand 
.sized  mica  flakes  were  added  to  a  small  quantity  of  filtered  distilled 
water  to  ol)tain  a  slurry.  A  drop  of  slurry  was  picked  up  on  a  fine  brush 
and  the  drop  barely  touched  to  the  top  of  tlic  miniscus  of  the  distilled 
filtered  water  in  the  cuvette,  thereby  introducing  mica  flakes  into  the 
measuring  app.iratus.  Other  methods  of  sample  introduction,  including  use 
of  an  eye  dropper,  were  tried,  but  tliey  tended  to  set  up  convection  within 
tlie  cuvette. 


Kesulting  settling  velocities  wort'  then  comp.ared  with  the 


thcorct ical 


settling  rates  generated  from  the  equ.itions  of  Lerman  and  others  (1974) 


Komar  and  Keimer.s  (1978).  In  tlie  c.i leu  l.i  l  i ons  of 
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for  ;:iic.i  ol  2.9  ;',/cm  was  usetl,  midw.iv  in  the  nor 


settling  rates  a  density 
al  range  for  hiotito  and 
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1  L  .  ily.niTiinat  ion  oi  the  mica  by  standard  optical  techniques  showed 
that  till'  :;ica  traction  in  tlie  cores  was  composed  of  muscovite  and  a  lesser 
af.oun.L  ai  hiotite.  Finally,  we  compared  tlie  measured  settling  rates  of 
tsi-  :.ica  sitii  tlioso  for  quartz  spheres. 


SILVPE  AS  A  FACTOR  OF  GRjMN  SETTLIN'G 
hiiape  has  been  recognized  as  an  important  factor  in  the  analysis  of 
,-n.'U t:'_eii t ev  iiydraulic  methods  (principally  settling  tube)  for  over  100 
vi-ars.  viibbs  and  o  tilers  (  19  7 1 )  ,  Lerman  and  others  (1974),  Komar  and 
ileimers  (l'*7S)  and  Brezina  (1979)  have  summarized  the  development  of 
tatak in.;  concerning  the  hydraulic  importance  of  shape  and  have  contributed 
to  tl:e  formulation  of  equations  for  settling  velocity  which  take  grain 
si’.ape  into  account.  In  a  series  of  recent  articles  Baba  and  Komar  (1981  a 
and  bl  itave  begun  to  e.xamine  sliape  offt.'cts  of  various  types  of  natural 
particles. 

Lerr.an  and  others  (1974)  modified  the  equations  for  the  settling  of 
a  disc  of  "no  thickness"  in  the  Stoke's  range,  developed  by  Lamb  (1932), 
Payne  and  Pell  (1960),  and  Brenner  (1964)  by  adding  a  term  for  disc 
tiiickness.  Tiie  resulting  formulae  for  the  two  major  orientations  of 


fall  .ire: 


U  = 


B(|>  -  p)  q,  r"  .  , 

s  h  (edgewise) 


3.396  M 


B(p  -  p)  q  r 
s  h 


(1) 


(2) 


U 


5.  1  p 


(broadside) 
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WiK'ro:  b'  =  soLtiin;’  sjxied 

r  =  di.sc  radius 

h  =  Chickness 

y  =  dynamic  viscosity  of  water 

=  particle  dens i tv 
p  =  density  of  water 

g  =  acceleration  due  to  gravity  (981  cm/sec”) 

Komar  and  Reimers  (1973)  approached  the  effect  of  shape  on  settling 
velocity  by  scale  modeling  wit'n  peblt.  os  in  a  glycerine  settling  medium, 
tne  results  Dcing  equivalent  to  quartz  sand  and  silt  in  water.  They 
found  experimentally  that  the  Corey  Siiape  Factor  (CST)  introduced  by  Corev 
(1949),  Malalka  (19-i9),  Nca<‘'wn  and  Malai-ta  (1930),  and  Mcllown  and  others 
(1951)  gave  the  best  prediction  of  shape  effects  for  the  pel'hles  they 
used  in  their  experiment.  Based  upon  their  exper it.ental  results,  Komar 
and  Reimers  (1978)  developed  an  empirical  formula  for  si'ttling  velocity 
in  the  Stoke's  range,  taking  into  account  shape  effoet; 
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■.iuTo : 


1) 

'I 

U. 


CSF 


h 

:  (CSI-) 

i (CSF ) 


2/3 

(liD.D^) 

intermediate  axi.s  particle  diameter 
principal  axis  particle  diameter 
ii/ 

1  1 

small  axis  particle  diameter 
-0.378 

0.946(CSF)  when  0.4  CSF  0.8 

2.18  -  2.09(CSF)  when  CSF  <  0.4 


ll.ev  tur tiler  found  that  they  could  extend  the  range  of  the  grain 
.'.ices  !\)ilewing  Stake's  settling  to  Reynolds  numbers  of  up  to  0.10  or 
(or  particles  of  up  to  100  um,  and  they  were  able  to  empirically  extend 
il'.eir  data  to  cover  grain  sizes  up  through  pebbles. 


Fo 

is 


r  particles  like  mica  where  h  is  small  relative  to  D.  and  D 

1 

0.4  and  the  general  equation  (3)  becomes: 


r 


CSF 


U  -  (1/18  u)  - - (p  -  p)  gdiD.D 

[2.  18  -  2.09  (h/.^D.Dj)  ]  ^  ^ 


Ihe  particie  thickiu'ss  h  may  bt.’  expressed  as  a  ratio  with  respect 


to  U,  i.i'.,  D/100,  i)/50,  D/25  ....  ])/ ^  and  therefore  in  the  ranges  in 


wiilch  are  workinp,,  the  wquations  of  Lerman  and  others  (197  <)  and 


e(iuation  (3)  are  similar.  Fur  example,  let  D.  =  D.  =  D,  x  =  ^  =  and 

‘  h  h 


r 


1 
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riion  otiu.it  ion  (4)  becomes: 


U 


_1_ 

18n 


1 


(2.  IS  -  2.09  -) 
X 


(P 


u  = 


(2.  tS 


_ 

2.09  -) 
X 


Stokes  Forr-.ula 


Our  r.onsuremonts  of  It  for  mica  flakes  yielded  a  ranye  of  h  0.007D 

1 


0.05  D.  For  x  =  1>/1>  = 


0.01  =  100,  and  subst itutin:4  in  (5),  we 


obtain: 


U  =  0.021 


ISii 


(p,  -  p)  gO" 
S 


Lerman'.s  formulae  are  of  the  general  form 

0  =  g -  p)  qj^r" 

cu 


wliere  c  =  3.396  to  5.1.  Sub.stituting  as  before  we  obtain 


U  =  g(p  ,  -  p)  0" 
2  c;ix 


Mu  1 1  i  p  ly  i  a;.;  bv  IS/ IS  yields 


i: 


^  -  p)  gif’) 


I'-lek.s  I'orr.uia 


-or 


'.ct  I'.in '  s  '  orriul.i  to  be  equivalent  to  Komar  and  Keimers  at  x  “  100 


•  e  -I't 


J_8 

') 

•-  >■  A 


0.  021 


18 _ 

2(100)  (0.021) 


c  =  A . 28 


ihis  value  lies  befaeen  Lerman's  tvo  values  of  3.396  and  5.1  and  therefore 
•alaes  uerived  from  the  tw('  sets  of  formulae  will  be  close  to  each  other 
ivi'r  the  ranees  we  are  examining. 


AssLL-nias;  D  =  =  D  (where  D  is  an  average  diameter)  is  an  idealized 


car.e  wiiica  all  mica  flakes  do  not  meet.  Anv  deviation  from  D  =  D,  will 

1  1 


oiivloualy  cause  some  deviation  in  the  calculation  of  y , 
parameter  >  =  0/h  becomes  the  shape  parameter  directly  related  to  the  Corey 


Mien  D.  =  D, ,  the 
1  1 


.'diape  f.ictor. 
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TUI-  HYDRAULIC  LQUIVALEXCV  OR  MICA 

In  our  oxpor imcnCs  micu  flakoa  seldom  settled  either  exactly  broad¬ 
side  or  edyewise.  They  usually  assumed  an  orientation  between  the  two 
extremes,  neither  perpendicular  nor  parallel  to  the  gravitational  vector 
whicii  t)iey  maintained  tlirougliout .  Stringham,  et  al  (  1969)  also  showed  that 
particles  settling  at  low  Revnolds  numbers  maintain  their  initial  orien¬ 
tation  which  need  not  be  perpendicular  to  the  settling  direction. 

Figure  2  shows  Llie  experimentally  determined  settling  velocities  of 
mica  flakes  plotted  on  the  family  of  curves  of  Lerman  and  others  (1974) 
and  Komar  and  Reimers  (1978)  for  thicknesses  of  U/2C0,  D/50,  D/20,  and  D/10. 
Komar  and  Reimers  (  1978)  equation  begins  to  diverge  from  those  of  Lermar. 
(1974)  only  at  thicknesses  of  D/20  and  greater.  The  thickness/D  ratios 
of  32  mica  grains  from  a  split  cf  the  slope  sediment  used  for  the  settling 
velocity  experiicent  were  determined  by  scattering  some  quartz  beach  sand 
on  an  Sill  stub  and  then  sprinkling  the  mica-rich  sloi'e  sample  over  it. 

Some  of  the  mica  flakes  landed  on  edge  and  we  were  thus  able  to  measure 
thickness  to  D  ratios.  Thicknesses  ranged  from  0.7%  of  D  to  5%  of  D, 
with  tile  average  at  2%;  and  a  standard  deviation  of  1.15%.  Figure  2 
also  siiows  tiiat  most  of  tlie  velocities  of  mica  flakes  wliicii  we  measured 
fall  within  tiie  envelope  created  by  the  theoretical  curves.  Scatter  at 
tlie  coarser  eiul  is  probalily  mostly  duo  to  variation  in  thicknoss,  scmin-.otrv 
(c.!;.,  D^  D^)  and  in  settling  orientation,  duo  to  tho  fact  tiiat  some 
particles  in  tiiis  size  raugi-  are  beginning  to  exceed  tiie  Stake's  low 
Reynolds  numijor  constraint  on  iIk-  Stukos  equation. 
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.\L  .-oai'scr  (.'iid  of  tlio  sea]  o,  par  t  ic  los  fall  closer  to  Lernan  and 

oLiier  .  U  '74)  equations  than  to  Komar  and  Reimers  (  1978).  Figure  3  shows 
tile  ui.ii'.eter  of  spherical  quartz  particles  of  density  of  2.65  which  would 

.i'ti!e  u  rates  of  the  variously  sized  mica  flakes  determined  from 
4et::an  inu  .)Lhers  (  1974)  equations  for  h  =  D/200,  D/50  and  D/20.  Our 
e.\;'.e  r ;  r.i-iu  s  sliow  tliat  mica  ranging  in  diameter  from  62  ym  to  250  ym, 

•  iir.  a  tiiickness  of  no  more  than  5%  of  the  sieve  diameter,  is  the 
r.vdraui  ic  L  quivalent  of  5  ym  to  82  ym  quartz  spheres.  Very  fine,  sand 
sir.i'd  -'.lea,  tiien,  is  tlie  hydraulic  equivalent  of  silt  sized  quartz 
-p.'ieres  '..iiile  fine  sand  sized  mica  is  the  equivalent  of  silt  to  very 
fine  s.iiivi  sized  splieres.  Coarse  silt  size  mica  flakes  are  tlie  equiva¬ 
lent  o:  fine, and  very  fine  quartz  silts. 


C0.\CLUS10.\’S 

1.  Ihe  equations  for  determining  settling  velocity  of  a  disc  pre- 
.-ienti'd  l)y  I.erman  and  others  (1974)  and  that  developed  by  Komar  and 
iieimers  (  1978)  for  particles  wliere  the  Corey  Shape  Factor  <  0.4  are 
.similar  lor  thicknesses  of  D/50  or  less. 

2.  hzper iments  on  fine  sand  to  coarse  silt  sized  mica  flakes  using 
a  ho  I  o  .’.r.ipn  ic  micro-velocimeter  showed  tliat  grains  mo.st  often  do  not 
orient  themselves  perpendicular  to  the  flow  field  as  they  settle  but  that 
thev  tend  to  settle  at  some  orientation  between  broadside  and  edgewise. 

3.  .'lica  in  the  coarse  silt,  very  fine  sand,  and  fine  sand  sizes  is 
the  hvdr.iulic  equiv.iU'iU  of  silt  .ind  vi-rv  fine  sand  sized  tjuartz  spheres 
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acco  i\i  i  r. lo  l1k‘  reJ  at  iunsliips  shown  in  Figures  2  and  3. 

Contrary  to  provions  assumptions,  mica  in  the  coarse  silt  to  fine 
saiui  siaos  is  not  the  hydraulic  equivalent  of  clay. 
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FIGURE  CAPTIONS 


Holographic  micro-veiocimeter  set  up  used  to 
measure  the  settling  velocity  of  the  mica  flakes. 

Figure  2.  Lcrman  and  others  (1974)  equations  for  the  settling 

of  discs  both  broadside  and  edgewise  plotted  along 
with  that  of  Komar  and  Reimers  (1978)  for  thickness 
of  D/200,  D/50,  D/20,  D/10.  Note  that  the  equations 
yield  values  which  are  very  close  to  each  other  at 
thicknesses  of  D/50  or  less,  thicknesses  most  prev¬ 
alent  for  mica  flakes.  X's  are  values  of  mica  flakes 
we  measured. 


3.  The  hydraulic  equivalency  of  mica  to  quartz  spheres. 

Mica  sizes  are  coarse  silt  to  fine  sand  and  thick¬ 
nesses  are  D/200,  D/50,  D/20. 


Figure 
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In  Situ  Holographic  Measurements  of  the  Sizes 
and  Settling  Rates  of  Oceanic  Particulates 

Ki  ndafi.  L.  Cardfr.'  RoBi  RrCj.  StfAVARi).  and  I’ftir  R.  Be  tzi  r 

Department  al  Marine  Setenee,  i  nteersite  at  Siitith  hlanja.  St.  Feter^htira.  tlarija  7.00/ 


A  froe-llo.umg  scilimenl  trap  equipped  »ilh  a  holographic  particle  tclocimeier  iHI'Vi  was  dcpioted  for 
14  4  hours  al  a  depth  of  .M)  ni  in  tlie  western  North  Atlantic  Ocean  The  ststeni  recorded  the  in  situ  sizes, 
shapes,  orientations,  and  seltline  rates  of  microscopic  particles  moving  thriiugh  the  laser  beam  The 
priniarv  data  reduction  revealed  panicles  from  the  system's  lower  limit  ol  resolution.  15  micrometers  in 
diameter,  to  250  micrometers  in  diameter  with  settling  velocities  ranging  from  OUlW  to  1)  2.502  cm  s 
1 16  lus  Ill  day  I  Individual  panicle  densities,  calculated  from  a  moditied  Stokes  equation,  ranged  Irom  I  .'7 
III  5  It)  g  ml.  The  presence  lU'  high  density  particles  was  independciilly  corroborated  through  indivijual 
particie  .inalvsis  ol  the  trapped  niaieri.il  with  a  computer-controlled,  scanning  electron  microicope 
equipped  with  an  enerey  ihspersive  -\-r.iv  .inaly/er  In  the  luiure.  in  situ  holographic  systems  might  i’.i.  used 
to  further  our  undersi.mdine  ol  primarv  productivitv.  sediment  erosion  deposition,  and  p.triicle  aggrega¬ 
tion  disruption  diss-  *unon 


ImRIIDI  l  IKIN 

Si/e.  shape,  and  density  are  key  variables  controlling  the 
settling  dynamics  and  late  of  ni.irme  particulate  matter.  The 
classical  method  of  determining  si/e  and  shape  is  to  collect 
suspended  particulate  matter  with  .i  device  such  ts  a  Ntskm 
bottle  l(iar,/iin.  iy7l)|  or  sediment  trap  and  then  concentrate 
the  particles  onto  a  memhr.ine  niter  or  in  a  settling  chamber  for 
microscopic  analysis.  The  uniyersal  problem  with  this  ap¬ 
proach  Is  determining  the  degree  to  whieh  the  si/e  distribution 
and  shapes  have  been  altered  by  the  collection  and  con¬ 
centration  sequences,  /■•niniier  fl’l-S]  .md  /rent  el  al  [in7s) 
have  show  n  that  sample  handling  can  cause  particle  disriiptiim 
Hi\lt"p  anii  l.JnuntJ  (  |vr6j  .md  /fiv/iop  |  I'l"]  h.iye  used  a  high 
volume  in  situ  pump  to  isol.ite  microscopic  oceanic  suspended 
matcri.il.  hut  even  with  this  well-engineered  system  there  is  ,i 
possibility  ih.ii  friable  material  will  not  maintain  ns  integrity 
through  the  scp.iration  process  There  is  considerable  adsan- 
tage  in  using  p.issne  m  situ  systems  to  characieri/e  seillmg 
panicles,  as  the  .tbiue  uncertainties  are  obyiated 

Pi  rposi 

.A  holographic  particle  velocimeier  (MP\T  has  been  lievel- 
oped  in  our  l.iboratory  that  measures  si/e.  shape,  and  seillmg 
velocity  of  mdivulual  particles  [<\iri/ir.  I')'');  CarJer  aiul 
Meters.  I‘)7‘)]  The  l.ibor.itorv  version  of  this  device  has  been 
modilied  for  use  in  a  frec-tloaling  sediment  trap  for  in  situ 
studies  the  purpose  of  this  report  is  to  describe  the  holo¬ 
graphic  device  and  present  some  results  from  what  are  believed 
In  be  the  lirsi  in  situ  settling  me.isurements  ol  individual,  micro 
scopic.  and  oceanic  particles. 

riiniRV 

riic  technique  is  an  applic.ilion  of  far-lield  or  FrauiiliolVer 
holography  |  I  hanipsun  et  a/  .  Idb’]  The  in-hiie  hologram  is  a 
record  ol  ihe  interlcrence  bciween  the  lar-lield  dilfraclion  pat¬ 
terns  sc.iitered  by  the  p.irlicle  and  the  collinc.tr  background 

N.'w  on  ',.ifsr.iin..it  .It  ( ITi.c  ol  Sn.KC  .iiiJ  lcrrcstrt.il  AppiK.itioiis, 
N  AS  y  VV  .ishrncion.  1 )  i  2'’24r. 

Copyright  l'i82  bv  the  .\iiiciic.in  t icopin  su.il  I  iiioii 

IVipcr  niinibcr  2t  oi'i| 

III  Js.1122'  s2  illl2l  o  l'l|  ys  no 


[Carinnahl  et  al.,  IdsO:  Tlinmpsan  el  al..  1‘Jb'']  Each  holcsgr.im 
provides  a  permanent  record  of  the  three  dimensional  position 
and  geometric  cross  section  of  all  particles  in  the  laser- 
ill  ii  mi  ii.ited  sample  volume  I. '.27  mil  I  Figure  1 1. 

The  laboratory  HP\'  was  modilied  three  ways  for  use  at  sea. 
First,  the  spatial  filter  was  removed  since  its  critical  alignment 
might  not  h.ive  wiihsiood  the  stresses  associated  with  depiov- 
meiit  .ind  recovery  of  the  sediment  trap  array  Second,  the 
imaging  lens  was  changed  from  50  to  72  mm  focal  length. 
Fin.illv.  the  I  mW  HcNe  laser  w.is  upgraded  to  2  niAA  The  last 
two  modilications  were  eiVccted  to  opiimi/e  the  device  for  faster 
settling  particles.  Fliat  is.  with  reduced  magnilic.ition  a  larger 
sample  held  was  recorded  .ind  with  incre.ised  laser  intensity  a 
shutter  speed  f.ist  enough  to  stop  particle  motion  on  lilm  was 
achieved.  To  d.impen  all  motion  other  than  that  due  to  gravi¬ 
tational  settling,  the  sample  volume  was  isolated  in  j  .-ii-cin  tall 
settling  chamber  I  Figure  1 1.  The  4  5  >  4.5  cm  square  chamber 
contained  two  sets  of  how  dampers,  one  5  cm  long  and  the 
other  III  cm  l-.ich  of  the  sets  consisted  of  about  50  luvtaposed 
cylinders,  each  i)  .5  cm  in  di.imeter. 

The  111  situ  UPS'  was  remotely  controlled  with  a  digital  timer 
progr. lined  to  trigger  evposurcs  at  '.5.  4.0.  and  15.5  s  after 
starling  the  laser  The  lime  between  evposurcs.  the  number  of 
evposurcs  sequence,  and  the  number  of  sequences  hour  were 
programabic  variables.  With  the  present  lilm  pack,  the  system 
can  collect  up  to  2511  cvpiisurcs. 

E\l‘l  RIMI  M 

The  Hl’\ -equipped  sediment  trap  collected  for  14  4  hours 
starling  from  I'HiO  June  I'l.  I'lSO.  at  a  site  i2b  N.  'ib  W’l  .iboul 
fill  nautical  miles  east  of  the  H.ihamas  I'lie  components  de¬ 
picted  in  Figure  1  were  mounted  in  the  ni.iin  cylinder  |4|icm  in 
dianieleri  of  the  Ir.ip  itself  The  ir.ip  w.is  located  at  'll  m  depth 
on  a  Iree-lloaiing  .irr.iv  composed  ol  four  parts  connected  by 
polyeihvlene-slic.iihed.  steel  c.ible  iFigure  2l  In  .iscending 
oriler.  these  were  ill  the  sediment  trap  with  si, unless  steel 
protective  c.ige,  |2|  ,i  set  of  three  subsuiMce  ho. its  starting  5  m 
above  Ihe  trap,  i.'-i  .i  set  of  three  surl.ice  boats,  and  i4i  .i  suriace 
m.irker  buoy  equipped  with  .i  strobe  light  .ind  r.idar  rcllcctor 
T"  d.impen  vv.ive-indiiccd  vertical  motion,  muiiiple  stranqs  of 
shock  cord  1 1  II  cm  diametcri  were  strung  between  each  of  the 

llo.lls 

Nf)\  \  scuba  divers  were  unable  to  detect  any  motion  or 
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I  lu  1  \n  im-.L-alcd  Nciiemaiic  iliacram  of  iht  in  >ilu  holojiraphic 
parii^ic  M'locimclLT  (HI*V  i  components  that  were  housed  in  the  sedi- 
nieni  trap  The  square  settline  lulv  was  open  tiniy  at  the  top  and  had 
two  seis  of  Mow  dampers  to  mtmmt/e  water  motion  in  the  sampling 
path 


harmonic  iibratiotts  in  the  arrai  just  prtor  to  tis  recoiery 
Neither  was  ans  trap-induced  lerttcal  motion  detected  with  the 
holographic  records.  That  is.  the  sanation  in  settling  velocilv 
per  particle  was  less  than  that  due  to  the  precision  ol'succesiie 
frame  registration  l  -  I  pml.  This  work  was  conducted  under 
moderate  sea  ctmditions  iwaie  heights  of  0  .a  I  m  and  winds 
less  than  HI  knolsl,  •Xddiltonal  testing  of  this  array  is  required 
to  learn  if  it  can  minimi/e  lerlical  motion  of  the  trap  under 
sesere  surface  conditions. 

Ml  IflllDS 

I'pon  recoiery  of  the  trap,  the  holographic  film  was  desel- 
oped,  ,ind  the  particulate  matter  from  the  bottom  chamber  of 
the  tr.ip  w.is  concentrated  onto  pre-tared,  0.45  jmi  pore  diam¬ 
eter.  membrane  Miters.  To  determine  the  total  mass 
tlus  m'  ilai,  the  tilters  were  weighed  at  our  shore-based  labora¬ 
tory  In  addition,  portions  of  the  pads  were  analyzed  by  a 
computer-controlled,  scanning  electron  microscope  equipped 
with  an  energy  dispersise  N-ray  analy/cr  at  the  State  l.'niver- 
sity  of  New  VUrk  SyractiseiD  L.  Johnsoni 

The  holographic  data  were  also  analyzed  at  our  shore-based 
facilities  n.ita  reduction  was  performed  manually  by  recon¬ 
structing  the  holographic  images  onto  a  white  screen  at  fItX) 
A  magnilic.ition  and  measuring  the  particle  dimensions  and 
displ.icemeiit  between  Iraitie  with  a  sernier  micrometer  1  >  5 
pinl  rite  density  ol  e.ich  particle  was  calculated  based  on  the 
settling  selocitv  and  loluine  estimated  from  the  two- 
dimension.il  image. 

The  accuraev  and  precision  of  the  data  reduction  techniques 
were  tested  m  two  laboratory  settling  esperiments  using  the 
s.inie  hologr.iphic  componeiits  as  deployed  in  the  lield  The 
esperimenis  were  conducted  with  presi/ed  (-t»  45  jiml  calcitc 
and  sphalerite  grains  Two  methods  ol  estimating  the  p.irticle 
\ Illume  were  compared  I  he  lirst  used  the  aserage  of  seser.il 
measured  dimensions  from  the  protected  cross-sccliottal  image 
as  the  spheric.il  equiialeni  diameter  in  the  standard  Stokes 
settling  eqii.iiion  for  spheres  ( <7  <ii.  Id'4|  Mie  second 
Used  the  .ippropn.ile  me.isured  dimensions  oi  the  p.irticle  .is  the 
maiiu  and  minor  .ises  ol  ,i  spheroid  I  liese  axes  were  each  used 
Is  ihe  roi.iiional  .iMs  l.i  .ipproxim.ite  obi. He  .ind  prolate 
spheroids  Hi  .issuming  the  polar  .ims  w.is  either  parallel  to  or 


1 

perpendicular  to  the  grai national  axis,  two  densities  each  wei 
calculated  by  using  equations  in  Lcrman  cl  til.  [  1474]. 

Compared  with  the  published  values  for  each  mineral  [Hur 
bui.  1971],  the  spherical  approximation  gave  the  most  consis 
ent  density  :  2.49  -f;  (1.22  calculated  versus  2.72  g  ml  publisho 
for  calcite  and  .v66  r  10  calculated  versus  4,0x0  10  g  * 
published  for  sphalerite.  In  both  cases,  the  spherically  derive 
densities  were  within  10''..  of  the  published  values.  Pycnometo 
determinations  on  the  same  samples  gave  2  b5  x  0.22  g  ml  fat 
calcite  and  4  15  x  0  22  for  sphalerite.  Occasionally,  individua 
spheroidally  derived  densities  were  much  closer  to  the  pub 
lished  values,  but  the  variation  in  the  measurements  was  quit) 
high  las  much  as  .500'. ,l  reflecting  the  inability  to  measure  lh» 
third  dimension  with  a  single  la.ser  system.  The  spherically 
derived  densities  with  a  potential  error  of  about  10 "  seem 
adequate  for  this  initial  attempt  at  quantifying  these  properties 
in  situ. 

Ruui  TS 

The  primary  data  on  each  hologram  is  a  panicle  size  distri¬ 
bution  of  all  the  particles  in  the  sample  volume.  Figure  3  is 
ineluded  as  an  example  of  such  data.  The  30  particles  depicted 
in  this  ticurc  represent  all  particles  in  excess  of  15  urn  that  were 
recorded  and  reconstructed  from  a  single  hologram  .^t  present. 
It  IS  not  feasible  to  integrate  such  a  small  time  increment  to  get 
an  estimate  of  mass  llu.x.  Under  developement  is  a  semt- 
automated  data  reduction  system  that  will  facilitate  the  gener¬ 
ation  of  a  sufficient  amount  of  optical  data  to  permit  estimation 
of  the  nitiss  ilux.  Through  gravimetric  analysis  of  the  collected 
material,  the  mass  flux  at  30  m  was  estimated  to  be  149 
mg  m-  day 


1  ig  .'  \n  unss.ilcd  scnem.iiic  .li.ier.im  .'f  ihc  Iree-Mo.uing  sediment 
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The  terminalion  iif  ihe  si/e  disinbuiion  at  15  //m  diameter 
was  dictated  hy  the  resoiiitmn  of  the  in  situ  system  as  it  lacked  a 
spatial  filter  and  had  lower  maenification  SEM  1;D\A  analy¬ 
sis  of  the  trapped  material  revealed  that  ‘>3  ..  of  the  sample 
mass  was  contributed  by  particles  whose  diameter  cxceded  In 
;<m  \S  hile  there  would  he  obsious  advantaees  in  extendine  die 
resolution  of  the  system  to  smaller  si/es.  the  present  comieur- 
ation  accounted  for  the  bulk  of  particles  which  contributed 
siitniticanily  to  the  mass  tlux  at  this  oceanic  sue. 

Some  of  the  faster  settline  particles  were  selected  to  show  the 
ranee  of  si/es,  settline  selocities,  and  densities  found  diirine  this 
m  situ  experiment  i  f  able  1 1  The  velocities  raneed  from  dOl'fil 
to  I'T.'iH)  cm  s  or  dailv  excursions  from  16  to  I'b)  m.  Density 
calculations  U'r  these  same  materials  raneed  from  1..'*  to  5.10 
i:  ml  1'  leurcs  4ii  4c'  are  photoeraphs  of  the  reconstructed  holo¬ 
graphic  images  of  three  particles  from  Table  I  having  densities 
of  4  Id,  ’,40.  and  I  .'7  g  ml,  respectively  figure  4i/  was  from  a 
deployment  of  the  1 1  l’\'-sedinient  trap  to  100  m  near  the  same 
site  1  days  later.  This  thin  vvalled  sphere,  tentatively  labeled  an 
invertebrate  egg  case,  was  included  as  an  example  of  one  of  the 
most  unusual  particles  recorded  during  the  experiment. 

DiSI  f  SSKIN 

The  results  show  that  it  is  possible  to  capture  in  situ  the 
hoKigraphic  images  (if  microscopic  marine  particles  and  use 
them  to  measure  their  individual  si/es.  shapes,  settling  velo¬ 


cities.  and  orientations.  ,As  far  as  we  can  determine,  these  are 
the  first  in  situ  measurements  of  these  variables.  Although  there 
are  many  possible  improvements,  the  success  of  this  experiment 
should  spark  the  application  of  this  technique  to  a  variety  of 
fundamental  oceanographic  research  problems  including  pri¬ 
mary  productivitv,  sediment  erosion  deposition,  and  particle 
aggregation  disruption  dissolution.  \  distinct  advantage  of  the 
system  described  herein  is  the  ease  with  which  it  can  be  altered 
to  study  selected  particle  si/e  and  density  ranges.  Magnification 
can  be  controlled  by  the  imaging  lens,  and  the  exposure  fre¬ 
quency  can  he  programed  to  capture  cither  slowly  or  rapidly 
settling  panicles 

One  ot  the  more  interesting  aspects  of  this  study  was  the 
density  approximations  for  some  of  the  particles  noted  at  this 
location  Three  of  the  1 1  listed  in  Table  I  had  densities  in  excess 
of  4  It  g  ml.  considerably  higher  than  one  would  predict  from  a 
know  ledge  of  typical  oceanic  suspended  particulate  matter  that 
includes  earbonates.  elay.  amorphous  siliea.  quart/,  and  organ¬ 
ics.  ,\s  this  IS  a  new  technique,  independent  verilic.dion  of  these 
results  was  sought  by  using  a  scanning  electron  microscope 
equipped  with  an  energy  dispersive  .\-ray  analy/er  and  image 
an.ilvMs  system  This  measured  particle  si/es  and  collected 
X-r.iv  data  for  71  elements  on  a  particle-by -particle  basis  from 
the  trapped  material  on  a  portion  of  the  tiller  pad 

Of  the  .S77  particles  cx.imined.  .34  |4  were  titanium-rich.  48 
til  I  were  iron-rich,  and  another  ,87  ( 10  ■  .i  exhibited  an  .X-ray 


TAHl.T  1  Individual  Suspended  I’.irlicul.ile  Si/c  ind  Seitling  X  elocilv  from  a  Holoei.iphie 
I’article  Veloeimcier-Sediment  Trap  Deploved  Dunne  Ihe  Niglil  ol  June  14.  |4SII 
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f  I;;  4  I’hi'liiL’rapIn  of  Iho  rccon>lrucIcd  holocraphu'  mumcA  of  llircf  panicli's  liMeJ  in  I'abk-  I  and  one  laree  vpheneal 
pariicic  toik-fiod  durme  a  later  e\pcrimeni  at  the  vame  mic  Ml  pa^^^ek•^  ate  oriented  with  the  itratitaponal  \eetor  poimine 
Iov^ar^l  the  hotlom  of  itie  paee 


>peeiriim  loiallt  l{:rcator  than  dS  '  i  made  up  of  iron  Thus. 
T'  (b\  iiiimberl  of  the  partieles  were  diiminaied  b>  elements 
norm.ilh  assoeiated  with  hiuh  density  mtnerals.  leg.  rutile, 
ilrneniie.  and  hemalilel  l  urlhermore.  the  elemental  distri¬ 
butions  of  these  particles  were  disitnet  ’roni  the  melallie  strtie- 
ttiral  eomponenis  of  the  setlimeni  trap  tihe  stainless-stcel  eaue 
and  I’.tlsani/ed  ehainl  Mlhoueh.  we  cannot  completely  rule 
out  eonl.iminalion.  the  above  evulenee  su.e.itests  that  we  cap¬ 
tured  hieh  ilensity  materials  setllini;  ihrotieh  the  upper  l.iyers  of 
the  oee.in  from  sources  other  than  the  trap  or  ship  As  three 
eruptions  ol  \li  St  Helens  volcano  occurred  vvithiit  .1  month  of 
the  cruise  illie  l.isi  vvilhm  '  davs  of  the  s.implini:i.  we  c.inui't 
rule  oiii  ihe  possibililv  that  it  provided  some  eoh.in  inputs  ti> 
lliesL'  samples 

In  .iildiiion  Ihe  p.irlicle-bv  p.irlivle  d.il.i  were  examined  to 
see  vvheiher  the  vlimensioivs  ol  the  hieli  ilensitv  particles  noted 


hoU'uraphically  were  consistent  with  those  measured  bv  the 
SIM  I  DXA  system  Several  titanuim  and  iron-nch  panicles 
were  found  iiiieroseopicallv  to  have  dimensions  essentially  co¬ 
incident  I  within  2  iim  each  of  Ihe  maior  and  niini'r  axes)  with 
Ihe  hieh  density  particles  reported  in  fable  I  Thus,  not  onlv  do 
the  compositional  characteristics  that,  when  translated  to 
oxides  fi'iind  III  sediments  t  fit).  1  efit).  1  eOl  match  the  holo- 
uraphicallv  ilelerminevi  viensiiies.  but  the  actual  dimensions  vsf 
Ihe  hiL'h  density  particles  recorded  bv  ihe  laser  are  the  same  as 
sever.il  recorded  by  the  scannmi:  electron  microscope  t)n  the 
b.isis  ol  this  corroboralion  and  the  resulls  of  our  own  l.ibora- 
torv  c.ilibi aiions.  the  deiisitv  .ippioxiniations  made  with  the  in 
siiii  1 1 1’\  .ire  reason.iblv  .iccur.ite 

I  he  calibr. lilt'll  expenmenis  vlenu'iisiraled  Ihe  impi'rlance  to 
the  spheroid, il  seitime  ev|iiatu'n  t'l  me  pri'per  setiline  I'rienta- 
lion  and  t'l  the  prt'per  rt'iation  of  a  tvvo-dinicnsional  imace  to 
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caicuLuions  of  a  ihrcc-dimcnsional  volume.  The  fiPV  simul- 
lanciuisK  records  the  actual  seltlme  direcuon  and  orteniatuvn. 
plus  the  cross-scciional  ^hape  of  each  particle  m  the  sample 
tield.  However,  this  single  laser  system  does  not  allow  deliniium 
of  an  a\is  of  rotation  needed  to  calculate  particle  volumes 
(mass).  Despite  this  limitation,  the  spherical  equivalent  den¬ 
sities  determined  for  two  standard  materials  were  wiihm  lu  . 
of  accepted  values.  Thus,  u  should  provide  similar  accuracies 
for  nu'st  compact  particles  le.g.,  dimensions  varving  by  less 
than  a  factor  of  2). 

On  the  basis  of  shape,  si/e,  settling  velocity,  particle  !I 
(Table  1  and  l-  igure  4i )  has  been  tentativeiv  labeled  as  a  fecal 
peiiei.  The  estimated  densitv.  I..O  g  nil,  is  consistent  with  wet. 
compact  organic  matter.  Despite  its  low  density,  the  si/e  id  this 
p<irticle  was  large  enough  lo  give  it  the  highest  settling  velocity 
measured  id.2.H)  cm  si.  Such  a  speed  is  cli>se  to  the  lower  end  of 
the  range  reported  for  fecal  material  from  the  equatorial  .Atlan¬ 
tic  \  liishur  I't  ill .  Id*”)  On  the  basis  id  this  estimated  rale  id 
descent,  the  particle  wiiuld  have  moved  though  the  euphoiic 
?one  in  less  than  a  dav.  This  underscores  the  potential  that  in 
situ  svsiems  have  in  addressing  biological  removal  processes 
that  mav  exert  primar>  Ci>nirol  i>ver  mass  tUix  in  Kdh  shallow 
and  deep  ocean  waters  [Lernum  er  a/..  19'^*’:  Snuinhi. 

Wiche  ft  ul..  19'’6;  lii.'^hup  et  a/..  1977;  I  imlfr  itiul  Snutll.  I9''2. 
fitinuuinJ  Rimuin.  \  97is.  Kiiatwr  ei  ai.  I97<>.  .Vi/iraaVr,  Id"*!] 

Our  experience  vvith  the  holographic  svsiem  indicates  that  it 
IS  cspecialiv  well  suited  for  analysis  id'  inorganic  maieriais. 
Primanlv,  this  results  from  their  being  more  refractive  ilian 
similarlv  si^ed  organic  material.  .As  a  result,  inorganic  material 
such  as  quart/  grains  and  clay  particles  have  high  Ci»nirasi 
dilTraciion  patterns  that  provide  a  sharper  hoK'grahic  image. 

The  Litilit  v  of  this  and  similar  devices  will  be  greailv  extended 
with  the  addition  of  a  rugged  spatial  filler  collimator  that  will 
improve  the  contra.si  of  the  recorded  dilfractii>n  patterns.  In 
iabi'ratorv  tests,  such  an  i>piical  arrangement  with  higher  mag- 
nificaih'n  and  spatial  tiltcring  has  alli»wcd  us  (i>  rcci>rd  and 
reciinsiruct  mineral  grains  as  small  as  2  /mi  in  diameter.  Such  a 
change  to  the  svsiem  will  also  be  an  advantage  in  studies  id' 
organic-rich  material  Mince  these  are  characterized  by  lower 
Cimtrast  ditVraction  patterns  (n  addiiii>n.  the  manual  image 
reduction  technique  could  he  automated  with  a  computer- 
based  image  anahsis  and  precision  frame  registraiiim  svsiem 
The  great  advantage  here  is  the  speed  and  accuracy  with  which 
the  data  from  large  numbers  of  holograms  could  be  processed 
W  ith  these  and  other  chanues,  passive  in  situ  optical  sxsiems 
may  help  increase  our  understanding  id  a  number  of  fundamen¬ 
tal  oceanic  processes. 

A(  kn<>\\  Ifilamcfifs  The  captain  and  crew  ivl  the  NO.A.A  vessel 
R  \  Kt'svnn  iwr  helped  make  (his  vvoik  possible  as  diil  (»coree 
Harvev  of  the  .Atlantic  Occani>graphic  and  Meteori'logical  l.ahoia 


lory.  We  especiallv  thank  Donald  Atwood  for  his  kind  invitation  to 
panicipalc  on  this  cruise,  This  work  was  supported  under  Ollice  of 
Naval  Research  contract  .N(MK)l4-75-(  -().''!^9. 
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.AhMract  Nuiricni  pro^lil.•^  I'rom  ihc  a>niincnial  shelf  of  the  northeastern  Gulf  of  Mexico  indicated 
considcrahie  near  hoiiom  enrichment  in  silica  and  nuraic  above  coarse  sediments  east  ot  the 
MissisMppi  della  In  contrast,  near  bottom  waters  of  the  carbonate  rich  West  Florida  Shelf  showed 
no  suen  enrichments.  Storm  related  suspension  apparently  produced  the  enrichments  because,  in 
near  Kutom  waters  south  I'f  Mobile  Bav.  silica,  nitrate  plus  nitrite,  and  suspended  load  increased 
suhNtatuiallv  as  a  winter  storm  front  passed.  Also,  laboratorv  simulation  of  resuspcnsion  by  stirnne 
(he  supernatant  Neawater  «>ver  a  clav  rich  core  produced  similar  increases  tn  silica  and  nitrate  plus 
nitrite,  with  amn  .a  being  the  apparent  precursor  to  the  nitrate  and  nitrite.  Most  of  the  nuineni 
increase  appeared  to  come  from  previously  deposited  sediments  in  the  early  stages  of  resuspension. 
L'sine  the  raiu>s  of  nutrients  released  to  sediments  resu.spendcd.  calculations  indicate  that 
resuspensum  of  as  little  as  I  mm  of  shelf  sediment  could  intcrmiiienliy  augment  overiying 
priKluctivitN  by  as  much  as  ItH)  to  200*’.*.  Thus,  resuspension  may  accelerate  nutrient  rccschnii  on 
continental  marems. 


iNTRODV  CT^ON 

As  i’\Kr  of  a  survey  of  the  continental  shelf  of  the  eastern  Gulf  of  Me.xico  during  midMay 
1974.  stations  were  occupied  m  two  adjacent  regions,  the  northeastern  gulf  shelf  from  the 
Mississippi  delta  to  Cape  San  Bias  (29°42'N.  85°24'\V)  and  the  West  Florida  Shelf  from 
Cape  San  Bias  south  to  Tampa  Bay.  At  each  station  (Table  I).  water  depth  was  70  m  or  less, 
and  nutrients  and  salinities  were  determined  for  surface,  mid  depth,  and  near-bottom  waters. 
The  resultant  protllcs  indicated  the  two  regions  to  be  different  in  an  important  respect  (Figs 
la.  b).  Near  bottom  waters  of  the  northeastern  gulf  shelf  were  considerably  richer  in  silica  and 
nitrate  than  those  of  the  West  Flonda  Shelf  even  though  the  near-bottom  salinity  profiles  for 
the  two  regions  were  virtually  identical.  The  empty  dashed  zones  in  Fig.  lb  show  where 
equivalent  enrichments  would  plot  if  they  also  occurred  in  the  West  Florida-Shelf  profiles.  But 
enrichments  were  only  found  on  the  northe.astern  gulf  shelf,  and  tho.se  for  stations  near  the 
Mississippi  delta  tended  to  be  the  highest  for  either  or  both  silica  and  nitrate  (Table  1,  Fig.  la). 
There  were  also  slight  increases  in  near  bottom  phosphate  on  the  northeastern  gulf  shelf 
compared  to  the  West  Florida  Shelf. 

The  enrichments  on  the  northeastern  gulf  shelf  might  have  resulted  from  four  proees.ses;  (a) 
transport  of  nutrient  rich  river  water  from  the  southern  United  States,  (b)  upwelling  of  deeper 
ofTshore  svaiers.  (e)  more  intense  near  bottom  regeneration  of  inorganic  nutrients  on  the 
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Table  1 .  Slalion  locations  for  survey  of  the  eastern  continental  shelf  of  the 
Gulf  of  Mexico,  15-18  ALav.  1474 


Station 

Latitude  ( N ) 

Longitude  (W) 

Shelf  section 

10 

29‘>36 

87°25 

NE  gulf  shelf 

11 

29‘>4I 

87°39.5 

NE  gulf  shelf 

12 

29<>46 

87°54 

NE  gulf  shelf 

13* 

29»57 

88°  14 

NE  gulf  shelf 

14* 

29°56.5 

88°23.5 

NE  gulf  shelf 

15* 

88°33 

NE  Euif  shelf 

A5’ 

29«I7 

88°26 

NE  gulf  shelf 

Ah* 

:9®20 

88°45 

NE  gulf  shelf 

C3 

29°  5  7 

87°  10 

NE  ftulf  shelf 

C4* 

29°33 

88°  13 

NE  gulf  shelf 

M7 

:9°4.t 

86°01 

NE  gulf  shelf 

M8 

29«44 

85°I4 

NE  eulf  shelf 

M4 

2905: 

86°15.5 

NE  gulf  shelf 

C2 

29°28 

85°49.5 

NE  gulf  shelf 

Ml 

27“45 

83°28 

West  Florida  Shelf 

m: 

27°52 

83°34 

West  Elonda  Shelf 

M3 

27<>5h 

83°43 

West  Florida  Shelf 

M4 

28°  21 

84°24 

West  Florida  Shelf 

M5 

28°29 

84°2I 

West  Flonda  Shelf 

Mb 

28°43 

84°20 

West  Flonda  Shelf 

Cl 

29°  1 .7 

84°02.5 

West  Flonda  Shelf 

•  Close  to  Mississippi  delta. 


SALINITY  SILICA  (/zM)  NITRATE  (^M) 


Fig.  1.  Caption  on  p  455. 
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northeastern  gulf  shelf,  or  (d)  an  injection  from  bottom  sediments.  Transport  of  river  water 
was  probably  not  responsible  because,  at  around  36  x  10~’,  the  near-bottom  salinities  were 
too  high.  Upwelling  was  also  unlikely.  Nutrient-salinity  plots  for  an  offshore  station  at  28°N. 
87°VV  (G.  Berberian  and  D.  Atwood,  personal  communication)  indicated  that  the  near¬ 
bottom  waters  at  Stas  1 1.  13.  and  15  had  much  higher  silica  concentrations  (13  to  16  pM) 
than  could  be  produced  by  the  upwelling  and  mixing  of  deeper  water.  Also,  if  upwelling  had 
provided  most  of  the  nutrients,  then  enrichments  in  silica  and  nitrate  should  have  occurred  at 
all  stations  having  upwelled  water.  But  near-bottom  waters  at  Stas  1 1  to  13  were  enriched  in 
silica  but  not  in  nitrate,  while  near-bottom  waters  at  Stas  14  and  15  were  enriched  in  silica 
and  nitrate  (Fig.  la).  Thus,  the  pattern  of  near-bottom  enrichment  was  not  consistent  with 
upwelling  as  the  major  source  of  increased  nutrient  concentrations.  There  is  no  reason  to 
expect  that  near-bottom  nutrient  remineralization  on  the  northeastern  gulf  shelf  was  more 
intense  than  on  the  West  Florida  Shelf,  so  that  means  of  enrichment  was  also  improbable. 

Thus,  much  of  the  near-bottom  enrichment  on  the  northeastern  gulf  shelf  appears  to  have 
been  caused  by  sediment-water  interaction,  with  turbulent  resuspensior  of  shelf  sediments 
being  a  reasonable  possibility.  Because  such  a  process  could  be  important  to  nutrient 
recycling  on  continental  margins,  we  conducted  two  related  studies:  (1)  observation  of 
nutrient  increase  and  sediment  resuspension  during  the  passage  of  a  storm  front  on  the 
northeastern  gulf  shelf  and  (2)  measurement  of  releases  of  silica,  phosphate,  and  nitrogenous 
nutrients  during  a  laboratory  simulation  of  sediment  resuspension. 

METHODS 

Salinity  vvas  determined  with  a  Beckman  salinometer  or  a  Guildline  Autosal 
conductometric  salinometer  standardized  against  standard  seawater.  Dissolved  ammonia  was 
determined  automatically  with  a  Technicon  AutoAnalyzer  11*  using  the  method  of 
Grasshoef  and  Johannsen  ( 1972).  Silica,  phosphate,  nitrate,  and  nitrite  were  determined  by 
the  automatic  techniques  discussed  by  Fanning  and  Maynard  (1978).  The  beam  attenuation 
coefficient  due  to  particles  suspended  in  seawater.  Cp(m'').  was  measured  with  a 
Hydroproducts'  906-912  transmissometcr.  A  particle-free  value  of  is  zero.  Suspended 
particle  loads  were  determined  gravimetrically  on  Nuclepore  membranes  (Betzer.  Carder 
and  Ecgimann.  1974),  Oxygen  concentrations  were  measured  by  Winkler  titration. 

NUTRIENT  INJECTION  BY  STORMS 

Storm  generated  resuspension  and  nutrient  release  were  observed  SSE  of  Mobile  Bay  at 
Sta.  2639  (29°53'N,  88°12'W).  water  depth  30  m.  For  a  5-day  period  (20  to  25  February, 
1978).  the  station  was  occupied  at  anchor,  with  transmissometry  profiles  taken  every  2  h  and 
near  bottom  water  samples  taken  daily  for  determination  of  nutrients  and  suspended  loads. 
Water  samples  for  nutrients  were  immediately  filtered  and  frozen  for  later  analysis,  and  the 
filter  pads  were  weighed  to  get  suspended  loads.  Salinity  samples  were  taken  every  2  h. 

A  winter  storm  front  passed  by  between  21  and  22  February  (Fig.  2).  About  0100  on  21 


Fiji  1  Prul'iles  of  nuiricnts  and  salinity  from  ihe  northeastern  Gulf  of  Mexico  and  West  Florida 
shelies  in  May  IV74.  Station  numbers  from  Table  I  are  shown  next  to  the  bottom  most  samples  of 
the  sianons.  The  dashed  prollles  indicate  stations  where  msufllcient  sampling  at  intermediate  depths 
introduced  some  uncertainty  in  salinity  distributions.  The  dashed  reitions  arc  explained  in  the  text. 


956  Kent  a.  Fanning.  Kendall  L.  Carder  and  Peter  R.  Betzer 


Fia.  Vanaiion  m  physical  parameters  durine  passage  ol'  .1  NTS  u  inter  sinrm  ironl  at  Sta. 

SSF  111  Mohile  Bay  (2h'’5,''N.  SS‘’li'\Vi.  Tstp:  wind  speed,  wind  direction,  sea  state.  Bottom: 
scTMcal  prsstik  of  beam  auenuatton  coeiTicicnt  of  thetransmissomctcr  ic'p)  im  ' 


February  the  sea  state  was  and  the  wind  speed  began  to  shift  from  8  kn  out  of  the 
southwest  toward  much  higher  speeds  coming  more  out  of  the  north.  The  maximum  speed 
occurred  about  24  h  later.  23  kn  out  of  the  northwest  at  0100  on  22  February.  By  that  time 
the  sea  state  had  reached  4.  About  30  h  after  it  started  the  storm  passed,  the  wind  returned  to 
about  8  kn.  and  the  sea  state  dropped  to  3.  In  the  last  half  of  the  storm's  passage,  the  turbidity 
of  the  yvater  column  began  to  increase.  Near-bottom  values  of  the  beam  attenuation  coefficient 
of  the  transmissometer  (fp)  rose  from  1. 1  to  2.7  m  '  around  0300  on  22  February  and  tlnally 
to  more  than  3.3  m”'.  About  the  time  of  the  highest  wind  speed,  the  surface  value  of  r.  also 
reached  a  maximum. 

The  storm's  passage  was  also  reflected  in  the  relationship  between  nutrient  concentrations 
and  suspended  load  in  near-bottom  waters  (Fig.  3).  Both  suspended  load  and  concentrations 
of  silica  and  nitrate  plus  nitrite  decreased  between  20  and  21  February  (see  Discussion).  But 
the  onset  of  the  storm  was  folloyvcd  by  corresponding  increases  in  suspended  load  and  in 
concentrations  of  both  silica  and  nitrate  plus  nitrite.  Suspended  load  increased  7-fold,  and 
each  nutrient  concentration  more  than  doubled.  The  bottom  half  of  the  water  column  had 
salinities  of  34.5  to  35.5  <  10  '  on  20  and  21  February,  and  the  water  that  moved  in  at  the 
storm's  heiitht  on  22  February  was  slightly  more  saline  (35.5  to  36.1  <  10'')  as  well  as  being 
richer  m  suspended  matter  and  nutrients.  The  salinity  nange  is  within  the  range  to  be 
expected  along  the  northwest  Florida  continental  shelf  (see  salinity  proHles  for  the  upper  30  m 
in  Fig.  la). 
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Fij,  3.  Variation  of  silica  and  nitrate  plus  nitrite  with  suspended  load  in  near  bottom  waters  at  a 
station  SSF  of  Mobile  Bay  iSta.  2h3d.  ;9°53'N.  88°  12'\Vv  A  storm  front  passed  the  spot  between 
2 1  and  22  February.  1978.  Arrows  mark  the  time  eourse  of  the  changes. 


The  near  bottom  suspended  load  during  the  storm  was  higher  than  those  from  other 
seasons.  The  same  station  was  occupied  during  summer  (five  determinations  in  July  1976) 
and  fall  (five  determinations  in  November  1978)  by  Bktzer.  Peacock  and  Jolley  (1978). 
Their  samples  were  taken  5  m  above  bottom  with  a  301  Niskin  bottle,  and  their  averages  for 
summer  and  fall  were  170  (±103)  and  380  (±228)  pg  I'',  respectively,  both  of  which  are 
much  smaller  than  the  1850  pgl  '  during  the  storm  of  22  February  (Fig.  3).  Betzer  et  al. 
(1978)  also  made  five  determinations  of  near  bottom  suspended  load  at  the  station  in 
February  1978  independently  of  the  data  in  Fig.  3.  and  their  average  was  1070(  r9 10)  pg  1 
which  is  also  higher  than  the  fall  and  summer  averages.  The  standard  deviation  for  the  winter 
sampling  is  4  to  8  times  those  for  other  seasons.  Such  a  large  scatter  is  consistent  with  the 
suggestion  that  winter  storms  produce  increases  in  suspended  load  by  resuspension  and 
transport. 

The  coincidence  of  a  storm  passage  with  marked  increases  in  suspended  matter  and 
dissolved  nutrients  is  unmistakable  and  strongly  suggests  the  changes  were  related. 
Resuspension  could  produce  dissolved  silica  in  two  ways,  the  mixing  of  pore  waters  into 
bottom  water  and  the  solution  of  silicate  particles  during  and  after  resuspension.  Marine  pore 
waters  contain  much  higher  concentrations  of  silica  than  bottom  waters,  and  release  of  silica 
from  stirred  sediment  has  been  demonstrated  (Fanning  and  Sciiink.  1969;  Fanning  and 
PiLSON.  1971).  Nitrate-nitrite  enrichments  could  have  occurred  through  the  intermediate  of 
ammonia.  In  shallow-water  reducing  scdimcni.  interstitial  ammonia  concentrations  can  be 
several  millimolar  while  the  overlying  water  column  has  much  less  ammonia  (Matisoef. 
Brickfr  Holdren  and  Kalrk.  1975).  /Vmmonia  may  also  be  sorbed  to  deposited  sediment 
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{Rosfnfeld,  1979).  Resuspension  of  a  few  millimeters  of  reducing  sediment  could  cause  a 
temporary  increase  in  ammonia  in  a  shallow  water  column,  and  subsequent  oxidation  could 
then  produce  nitrate  and  nitrile. 

The  suspended  matter  associated  with  the  suspended- load  maximum  on  22  February 
(Fig.  3)  contained  a  considerable  fraction  of  clay — about  50%.  This  conclusion  was  reached 
from  its  aluminum  content  (6.3'’<))  and  its  Si:  Al  ratio  (3.4  :  1).  Such  a  low  Si :  A1  ratio  means 
that  much  of  the  suspended  silica  was  associated  with  aluminum  in  aluminosilicates.  The 
estimate  of  50'’(>  clay  was  made  using  a  mass  conversion  ratio  of  8  units  clay  to  1  unit 
aluminum,  based  on  mineralogic  determinations  on  the  suspended  matter.  Thus,  the 
suspended  matter  near  the  bottom  during  and  after  the  storm  not  only  had  greatly  enhanced 
concentration  but  also  had  a  composition  consistent  with  an  origin  in  a  clay -rich  substrate. 

Recent  work  indicates  that  the  bottom  under  Sta.  2639  has  many  sand-sized  particles 
(DOYt  i  and  Sparks.  1980).  However,  just  to  the  west  the  sand  content  decreases  markedly 
as  fine  grained  pro-delta  facies  of  the  Mississippi  River  are  encountered.  Thus,  the 
resuspension  and  nutrient  increases  observed  at  the  station  may  have  originated  in  the  west 
and  been  earned  to  Sta.  2639  as  the  storm  moved  eastward. 

STIKRI.sn  (  XPI  RIMI  s  r 

To  investigate  nutrient  release  by  resuspension  under  more  controlled  conditions,  we 
conducted  a  laboratory  resuspension  of  deposited  sediment.  Because  the  sediment  in  the 
water  column  at  Sta.  2639  contained  considerable  amounts  of  aluminosilicates,  the  deposit  to 
be  anitlcially  resuspended  had  to  be  clay-rich  and  could  not  come  from  beneath  Sta.  2639 
(see  above).  We  lacked  cores  from  the  Mississippi  delta,  but.  on  cruise  78-R-2  of  the  R.V. 
Kescarc/’i'r.  vve  were  able  to  obtain  clay-rich  sediment  iVom  the  continental  margin  north  of 
Venezuela.  Those  deposits  derive  from  the  Orinoco  river  and  are  typical  clay-hearing  deltaic 
muds.  Bv  gently  lowering  a  weight  stand  with  two  core  liners  attached  1 3  cm  apart,  two  cores 
were  taken  simultaneously  in  103  m  of  vvatcr  at  1  l°(X)'N.  62°0rW.  The  cores.  G.C.  17  and 
G.C.  18.  were  6.7  cm  in  diameter  and  were  sealed  in  nitrogen-tilled  plastic  bags  and 
refrigerated  to  retard  bacterial  activity  until  opened  and  processed. 

Core  G.C.  17  was  used  to  characterize  the  sediment.  At  room  temperature  (22°Ck  it  was 
sliced  in  a  nitrogen  atmosphere  and  the  slices  squeezed  to  remove  pore  waters,  which  had  high 
concentrations  of  interstitial  silica  |200  to  300  pM)  and  phosphate  (4  to  9  pMl.  Most  of  the 
sediment  vvas  reducing;  the  top  centimeter  had  15  pM  nitrate,  but  the  rest  had  much  less  than 
5  pM  nitrate.  The  sediment  was  brownish  on  lop  and  greenish  below  5  cm.  Results  for  G.C. 
17  were  assumed  to  apply  for  G.C.  18  as  well  because  the  two  were  collected  so  close 
together. 

Core  G.C.  18  was  used  for  the  resuspension  experiment  at  22°C.  approximately  the  iit  ^iiu 
temperature  ( 19°C).  The  old  supernatant  seawater  vvas  carefully  removed,  and  I  I  of  low- 
nutrient  surface  seawater  (salinity:  36.57  \  10  ’)  vvas  put  over  the  sediment.  The  seawater 
was  very  slowly  added  through  capillary  tubing  so  that  it  trickled  gently  down  the  liner  wall 
above  the  sediment  until  the  water  column  over  the  sediment  vvas  28  cm  high.  The 
emplacement  took  3  to  4  h.  The  seawater  was  allowed  to  stand  for  30  min  and  then  sampled. 
A  6  cm  wide  Tellon  paddle  vvas  placed  about  10  cm  above  the  scdimeni-waicr  interface  and. 
over  the  next  20  h.  vvas  rotated  to  resuspend  the  sediment.  The  rotation  was  slow  at  first,  with 
the  outer  edcc  of  the  paddle  having  a  t.ingcntial  velocity  of  9  cm  s  '.  At  intervals,  rotation 
speed  was  increased,  and  samples  of  the  suspension  were  taken  I'or  the  determination  ol' 
suspended  load  and  dissolved  components  i Table  2). 
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rahlv  1.  Stirwi^  speeds,  suspended  loads,  and  nutrient  concentrations  duhni^  the  resuspension  experiment  on 
Core  (J.C.  US  from  the  yenezueian  continental  shelf  62°Ol'\V).  Supernatant  seawater  was  initially 

nutrient  poor  with  36.57  x  10  ’  salinity 
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•  Times  \shcn  me  rotation  speed  was  increased. 


Erosion  of  ihe  sediment  was  slight  for  the  first  30  min.  Then  the  stirring  speed  was 
increased  to  33cms  '  (Table  2.  Fig.  4).  and  the  suspended  load  increased  rapidly  with 
stirring  rate  until,  between  190  and  200  min.  it  reached  20  g  I  '.  Stirring  was  continued  for 

another  16  li.  and  ultimately  about  40  g  of  sediment  were  eroded.  The  area  of  the 

sediment- water  interface  was  35.3  cm^  so  between  1  and  2  cm  of  sedimeni  had  been 
resuspended,  assuming  a  surface  sedimeni  porosity  of  60  to  S0''.i  and  a  dry  particle  density  of 
2.7  g  cm  '. 

Nutrient  release  seemed  to  be  in  three  stages  (Table  2.  Fig.  41.  but  the  first  stage  was  the 
most  une.tpected  and  perhaps  the  most  significant.  It  occurred  during  the  emplacement 
operation  before  any  stirring  was  begun.  The  seawater  trickling  down  the  core  liner  wall  had 
almost  no  nutrients  (Table  2);  yet.  after  emplacement  and  30  min  of  waning,  all  nutrient 
concentrations  had  increased.  The  ammonia  increase  was  from  undetectable  to  4.5  pM  and 
silica  increased  25  fold  Thus,  a  release  of  nutrients  is  shown  at  zero  stirring  time  in  Table  2 
and  Fig.  4.  w  ith  the  initial  release  of  nitrogenous  nutrients  being  nearly  half  of  the  final  release 
obsersed.  DilTusion  across  the  sediment  water  interface  did  not  provide  all  of  the  release 

observed  during  emplacement.  .Xgain  assuming  60  to  80'hi  porosity  in  the  sediment,  the 

increases  in  silica  and  phosphate  would  have  reuuired  all  the  interstitial  silica  and  phosphate 
in  the  upper  5  to  10  mm  of  sediment,  and.  during  the  3  to  4  h  emplacement  period,  diffusive 
exchange  could  not  have  occurred  over  more  than  I  to  2  mm.  Thus,  the  initial  release  came 
partly  from  a  small  amount  of  sediment  resuspended  during  the  emplacement.  No  release  was 
anticipated,  so  the  suspended  matter  in  the  seawater  was  not  determined  before  emplacement, 
and  how  much  of  the  21  mg  I  '  suspended  matter  before  stirring  (Table  2)  was  due  to  the 
emplacement  is  unknown.  Thus,  the  nutrient  release  per  gram  ol  emplacement  resuspended 
sedimeni  is  difficult  to  ascertain  .Assummu  that  the  emplacement  water  initially  had  no 
suspended  matter,  we  can  estimate  minimum  release  ratios  (Table  3). 
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The  sccdikI  staee  of  nutricni  release  eame  after  30  iiiin  of  stirnitu  when  the  stirnni;  speed 
was  inercased  lo  38  em  s  '  and  extensive  erosion  beitan  (Table  2.  f-ic.  4).  .Ml  nuirieiu 
eoneentr.iiions  inereased  rapidly  during  the  resuspension.  and  the  relative  inereases  in 
phosphate  and  siliea  were  higher  than  that  of  suspended  matter  (Tig.  4).  Ratios  ol  nutrient 
released  per  gram  of  eroded  sediment  vvere  ealeulated  for  the  period  of  35  lo  60  min  stirring 
lime  when  the  erosion  was  most  rapid  (Table  3).  The  release  ratios  arc  roughly  two  orders  of 
magnitude  less  than  those  for  the  lirsl  stage. 
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The  third  stage  of  nutrient  release  was  between  3  and  1 9.5  h  of  stirring  and  was  weaker 
than  the  other  two.  Silica  and  nitrogenous  nutrients  were  released  slightly  and  phosphate  may 
have  even  decrea.sed.  however  the  total  amount  of  resuspended  sediment  nearly  doubled 
(Fig.  4).  The  implication  is  that  most  of  the  releasable  nutrients  are  in  the  upper  few 
millimeters  of  sediment. 

During  the  three  stages,  release  of  nitrite  and  nitrate  appeared  to  be  related  to  that  of 
ammonia  (Table  21.  The  first  major  change  was  the  appearance  of  ammonia.  As  erosion 
commenced  and  continued,  ammonia  concentrations  remained  fairly  constant  and  then  rose 
to  about  7  pM  after  1  h  of  stirring.  Meanwhile,  nitrate  concentrations  increased  slightly  and 
then,  after  the  first  hour,  began  to  increase  as  the  ammonia  concentrations  began  to  fail.  After 
3  to  19  h  of  stirring,  nitrate  became  the  principal  nitrogen  compound.  The  pattern  of 
ammonia  release  followed  by  increasing  nitrate  concentration  resembles  the  classic 
description  of  the  regeneration  of  nutrients  from  phytoplankton  populations — 
ammonification  followed  by  nitrification  (Richards.  1965;  HaRVEY.  1955).  Oxygen 
concentrations  measured  at  3.2  and  19.5  h  stirring  time  were  4.0  and  4.4  ml  (STPll  '. 
respectively.  Thus,  ample  oxygen  was  present  to  oxidize  ammonia  and  bacterial  substrate  was 
ample. 

The  (lux  of  nitrogenous  compounds  during  stirring  can  be  compared  with  fluxes  estimated 
for  other  nearshore  sediments  (Table  4).  The  rate  used  for  the  experiment  was  that  from  when 
erosion  was  extensive  between  35  and  60  min.  The  flux  in  the  experiment  is  much  higher  than 
the  other  estimates,  which  are  either  for  long  time  periods  (Billen.  1978)  or  for  benthic 
chambers  that  prevent  resuspension  (Nixon.  Oviatt  and  Hale.  1976).  In  contrast,  a 
resuspension  (lux  associated  with  erosion  is  for  a  single,  dramatic,  and  fairly  infrequent  event. 


ruble  i.  ICsrimuU  d /luxes  uf  inoreanic  forms  of  combined  nUrogen  from  uirious  areas 
of  the  eominental  margin  under  di/fereni  enndilions 
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If  the  interstitial  concentrations  of  core  G.C.  1 7  represent  those  of  the  stirred  core,  the  pore 
waters  released  during  stirring  had  about  250  pM  SiO;  and  5  pM  POj  .  Assuming  a  sediment 
of  60  to  80'’o  porosity  and  a  dry  density  of  2.7  g  cm  \  the  amounts  of  interstitial  silica  and 
phosphate  available  were: 

8.3  X  10  to  2.2  >  10  '  g  SiO .  per  g  solid 
and 


2.6  X  10  to  7.0  «  10  g  PDi  per  g  solid. 
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But  the  data  in  Table  3  show  that,  during  rapid  erosion,  the  release  ratios  of  silica  and 
phosphate  were; 

2.3  X  I0~‘  g  SiO;  per  g  resuspended  solid 

and 

1.7  X  10'-  g  POj'  per  g  resuspended  solid. 

Thus,  during  the  second  stage,  the  interstitial  waters  appeared  to  contribute  only  a  small 
fraction  of  the  total  silica  and  phosphate  released.  The  implication  is  that  the  upper  few 
millimeters  of  continental-margin  sediment  may  contain  a  considerable  reservoir  of  nutrients 
attached  to  the  deposited  particles  and  releasable  by  resuspension. 

DISCUSSION 

The  information  presented  suggests  that  resuspension  of  continental  margin  sediment  is 
important  in  fertilizing  nearshore  waters.  Near-bottom  nutrient  enrichment  along  the 
northeastern  margin  of  the  Gulf  of  Mexico  can  be  explained  by  such  a  process  as  can  storm- 
related  changes  in  nutrients  and  suspended  load  south  of  Mobile  Bay.  Alabama.  Further,  a 
laboratory  study  of  deltaic  sediment  from  the  Orinoco  demonstrated  substantial  nutrient 
release  from  rather  gentle  disturbances. 

The  question  of  the  true  importance  of  resuspension  as  a  fertilizing  mechanism  arises. 
However,  we  wish  to  emphasize  that  resuspension  is  at  best  a  mechanism  for  recycling 
nutrients  between  water  column  and  sea  bed.  Any  nutrients  released  by  resuspension  could 
only  have  reached  the  sea  bed  by  the  sinking  of  some  form  of  organic  matter  from  the  water 
column.  The  ultimate  source  of  most  oceanic  nutrients  is  land  runoff.  However,  resuspension 
could  play  a  role  in  sustaining  high  productivity  on  some  continental  margins  by  accelerating 
nutrient  turnover. 

For  the  Gulf  of  Mexico,  considerations  of  the  importance  of  resuspension  to  nutrient 
recycling  involve  calculations  of  potential  production  from  riverine  nutrients,  storm  intensity 
and  the  distribution  of  fine-grained  sediment,  and  near-bottom  er  .chment. 

According  to  work  of  the  U.S.  Geological  Survey  (1980).  during  1977  to  1979  the 
Mississippi  River  contained  an  average  of  8.25  pg-at.  P  r‘  in  both  dissolved  and  particulate 
form  and  an  average  of  164  pg-at.  N  I  '  as  ammonia,  nitrite,  nitrate,  and  dissolved  and 
particulate  organic  nitrogen.  We  wili  consider  that  all  forms  of  nitrogen  and  phosphorus  can 
be  used  by  phytoplankton  and  that  all  of  the  river  runoff  to  the  gulf  from  the  United  States 
has  the  average  concentrations  obtained  for  the  Mississippi.  The  Mississippi  River  runolT  is 
1.54  X  I0''ld  '  (Hole.man,  1968)  and  the  runolT  from  other  United  States  rivers  is 
0.41  X  10'*  I  d  '  (El-Sayep  ei  aL  1972).  making  a  total  runoff  from  the  United  States  to  the 
gulf  of  2  x  10'’  I  d  '.  The  area  of  the  United  States  gulf  shelf  is  3.34  x  lO"  m*  (Manheim. 
1980).  Therefore,  if  we  assume  the  Redfield  ratios  apply  ( 106  C  :  16  N  :  I  P:  Richards,  1965), 
the  United  States  rivers  could  sustain  an  average  shelf  productivity  along  the  gulf  coast  of 
63  mg  C  m  ’  based  on  the  total  riverine  phosphorus  or  78mgCm  ’d  '  based  on  total 
riverine  nitrogen.  Yet  the  measured  productivities  along  the  United  States  gulf  coast  are  2  to  3 
times  higher:  100  to  200  mg  C  m  'd  '  (El-Sayed  ct  ai.  1972).  It  is  possible  that  the 
phytoplankton  of  the  United  States  gulf  coast  photosynthesizc  more  atoms  of  carbon  per 
atom  of  phosphorus  or  nitrogen  than  the  Redfield  ratios  predict,  but  the  higher  shelf 
productivities  may  come  from  recycled  nutrients.  The  possibility  appears  more  likely  because 
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all  of  the  riverine  nitrogen  and  phosphorus  used  in  the  above  estimate  may  not  be  usable  by 
phytoplankton:  sediment  resuspension  could  provide  the  additional  source  of  nutrients. 

An  estimate  suggests  the  potential  enhancement  of  primary  production  by  resuspension 
along  the  gulf  shelf  with  storms  playing  a  major  role.  If  1  mm  of  shelf  sediment  is 
resuspended,  and  the  sediment  had  80%  porosity  and  2.7  g  cm“^  dry  density,  its  resuspension 
would  introduce  540  g  of  solid  material  per  square  meter  of  water  column.  If  resuspension  of 
shelf  sediment  from  the  gulf  shelf  releases  approximately  the  same  quantities  of  nitrogen  and 
phosphorus  as  occurred  during  the  experiment,  we  can  multiply  the  540  g  m  *  by  the  release 
ratios  in  Table  3  and  the  Redfield  ratios  to  predict  the  potential  enhancement  in  production. 
The  appropriate  Redfield  ratios  are:  0.0795  mg  C  (pg-at.  N)“'  and  1.272  mg  C  (pg-at.  P)“‘. 
From  Table  3.  release  ratios  for  the  emplacement  stage  predict  enhancements  ranging  from 
more  than  6900  to  more  than  9000  mg  C  m'^:  the  ratios  for  the  extensive-erosion  staee 
predict  enhancement  of  40  to  120  mg  C  m  '.  So  resuspension  of  1  mm  of  sediment  could  raise 
primary  production  considerably,  especially  if  the  phytoplankton  fix  more  carbon  than  pre¬ 
dicted  by  the  Redfield  ratios  (Goldman.  McCarthy  and  Peavey,  1979).  Further,  there  are 
many  storms  to  resuspend  sediment  in  the  Gulf  of  Mexico.  Fine  grained  resuspendable  sedi¬ 
ment  covers  approximately  200.000  km"  of  shelf  from  Mobile  Bay  (Doyle  and  Sparks. 
1980)  to  Mexico  (Shidler.  1977).  and  much  of  it  is  under  50  m  of  water  or  less.  Intense 
storm  fronts  cross  the  area  from  northwest  to  southeast  an  average  of  four  times  a  month 
(Henry.  1979):  during  the  hurricane  season,  an  average  of  four  tropical  cyclones  and 
hurricanes  run  ashore  across  this  sediment  (Simpson  and  Lawrence.  1971).  Each  storm 
takes  1  to  2  days  to  pass  a  given  point  on  the  shelf,  implying  that  the  nutrients  could  be 
introduced  fast  enough  to  be  important  to  rates  of  productivity. 

There  are  published  references  to  the  importance  of  recycled  nitrogenous  compounds  from 
the  sea  bed.  Ho  and  Barrett  (1977)  concluded  that  ammonia  enrichments  in  near-bottom 
waters  of  Caminada  and  Barataria  bays  near  the  Mississippi  delta  in  March  and  May  of  1973 
were  produced  when  waves  agitated  sediments.  Ni.xon  (1981)  presented  budgets  that  show 
that  ammonia  emitted  from  benthic  communities  is  an  important  part  of  the  nutrient  balances 
in  coastal  waters. 

Our  work  on  resuspension  has  some  puzzling  features.  During  the  storm  passage  at  Sta. 
2639  south  of  Mobile  Bay  (Fig.  3).  suspended  sediment  increased  1620  pg  I"',  sibca  rose 
2.4  pM.  and  nitrate-nitrite  rose  4.6  pM  throughout  the  bottom  15  m  of  the  water  column. 
The  release  ratios  as  calculated  from  these  field  data  are  considerably  larger  than 
observed  in  the  stirring  experiment  (Table  3).  The  sediment  resuspended  during  the  first  stage 
of  the  stirring  experiment  could  only  be  estimated  as  a  maximum.  If  less  were  actually 
resuspended  during  the  emplacement,  the  release  ratios  for  the  first  stage  would  have  been 
higher,  more  in  agreement  with  ratios  from  Sta.  2639.  The  high  fraction  of  clay  minerals  in 
suspension  at  Sta.  2639  suggests  that  the  suspended  matter  may  have  originated  elsewhere, 
perhaps  to  the  west.  Dunne  transit  to  Sta.  2639.  some  of  the  suspended  matter  may  have 
settled,  but  any  nutrients  pre.sent  were  dissolved  and  may  not  have  vanished  from  the  water 
column  as  quickly.  The  ratio  of  released  nutrient  concentrations  to  resuspended  sediment 
loads  would  appear  to  have  increased  during  transit  and  would  have  been  higher  at  Sta.  2639 
than  immediately  after  a  resuspension  event,  either  in  siiu  or  in  the  laboratory.  Another  puz¬ 
zling  feature  of  our  work  is  the  data  from  20  February  (Fig.  3).  On  that  day.  the  water  at  Sta. 
2639  had  nutrient  concentrations  higher  than  those  during  the  storm's  passage  and  yet  much 
lower  suspended  loads.  Again,  the  explanation  may  be  the  sinking  of  resuspended  sediment 
away  from  the  soluble  nutrients.  In  winter,  storms  pass  that  area  every  2  to  4  days.  Water 
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may  have  been  carried  in  before  20  February  with  high  nutrient  concentrations  and 
suspended  loads  and  then  remained,  allowing  the  suspended  matter  to  settle  while  the  nutrient 
concentrations  remained  fairly  high.  However,  without  more  data  we  can  only  speculate 
about  events  prior  to  20  February. 


CONCLUSIONS 

Evidence  from  nutrient  profiles  on  the  northern  rim  of  the  Gulf  of  Mexico,  from  increases 
in  nutrients  and  suspended  matter  during  a  winter  storm,  and  from  laboratory  resuspension 
studies  on  fine  grained  shelf  sediment  strongly  suggest  nutrient  release  during  the  resuspen¬ 
sion  of  sediments. 

The  results  from  a  laboratory  release  experiment  on  stirred  sediment  can  be  used  in 
evaluating  the  potential  impact  of  resuspension  by  normalizing  the  nutrient  release  to  the  sedi¬ 
ment  release  during  stirring,  giving  a  release  ratio  (Table  3).  Application  of  those  ratios  show 
that  as  little  as  I  mm  of  resuspended  shelf  sediment  could,  with  subsequent  mixing  into  the 
euphotic  zone,  significantly  fertilize  the  overlying  waters.  The  elTect  is  an  enhancement  in 
nutrient  recycling  on  continental  shelves,  thus  providing  a  reasonable  explanation  of  why  shelf 
productivities  could  be  higher  than  predicted  from  fluxes  of  river  nutrients. 

The  potential  release  of  recycled  nutrients  from  resuspended  sediment  should  be  included  in 
studies  of  nutrient  balances  along  continental  margins. 
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